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Abstract. There is accumulating evidence that (fundamental) scalar fields may exist in Nature. 
The gravitational collapse of such a boson cloud would lead to a boson star (BS) as a new type of a 
' compact object. Similarly as for white dwarfs and neutron stars, there exists a limiting mass, below 

, which a BS is stable against complete gravitational collapse to a black hole. 

■ According to the form of the self-interaction of the basic constituents and the spacetime sym- 
' metry, we can distinguish mini-, axidilaton, soliton, charged, oscillating and rotating BSs. Their 

.| , compactness prevents a Newtonian approximation, however, modifications of general relativity, as 

^ ' in the case of Jordan-Brans-Dicke theory as a low energy limit of strings, would provide them with 

h— ^ I gravitational memory. 

I , In general, a BS is a compact, completely regular configuration with structured layers due to the 

anisotropy of scalar matter, an exponentially decreasing 'halo', a critical mass inversely proportional 
I I , to constituent mass, an efi'ective radius, and a large particle number. Due to the Heisenberg principle, 

. there exists a completely stable branch, and as a coherent state, it allows for rotating solutions with 

quantised angular momentum. 

Q , In this review, we concentrate on the fascinating possibilities of detecting the various subtypes 

$_H ' of (excited) BSs: Possible signals include gravitational redshift and (micro-) lensing, emission of 

^ ' gravitational waves, or, in the case of a giant BS, its dark matter contribution to the rotation curves 

. of galactic halos. 

^ ! PACS numbers: 03.75.F; 04.20; 04.25; 04.30; 04.40.D; 05.30.J; 95.30.S; 95.35; 97.10.R; 98.35.J; 98.52.N; 98.52.W; 98.58.H 
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I. INTRODUCTION 



In this review, wc will assume that fundamental scalar fields exist in Nature, that in the early stages of the 
universe — they would have formed absolutely stable soliton-type configurations kept together by their self-generated 
gravitational field. Theoretically, such configurations are known as boson stars (BSs). In some characteristics, they 
resemble neutron stars; in other aspects, they are different and, thereby, astronomers may have a chance to distinguish 
BS signals from other compact objects, like neutron stars or black holes (BHs). BSs can be regarded as descendants 
of self-gravitating photonic configurations called geons (gravitational electromagnetic units) and proposed in 1955 by 
Wheeler [307]. 

In this review, we shall mainly concentrate on the results which are not included in the first reviews from 1992 
[152,184,185] as well as later in the Marcel Grossman meetings in Jerusalem and Rome [212,215]. Since then, the 
number of publications investigating how BSs could possibly be detected have increased; we intend to focus strongly 
on these more recent research results. This could provide astronomers a handling on possible signals of BSs. All of 
these theoretical results are only a humble beginning of the understanding where and in which scenarios a BS could 
be detected. We hope other researchers may find this review stimulating for gaining new ideas or for refining older 
research. 

A second intention of this review is to distinguish clearly the different theoretical models underlying the label BS 
which could lead to different observational consequences. 

The first distinction is that the matter part of a BS can be described by either a complex or by a real scalar field. 
Then, there can be different interactions: (a) self-interactions described by scalar field potentials, (b) minimal coupling 
to gauge fields, the scalar field can be carrier of a charge (electric or hypercharge, e.g.). Moreover, the BS scalar field 
can couple, in standard general relativity (GR), minimally, or, in scalar- tensor (ST) or Jordan-Brans-Dicke (JBD) 
theory, non-minimally to gravity. In the latter case, if the strength of gravitational force is influenced by the JBD 
scalar field, there is an interaction of the BS scalar field with the real JBD scalar field; this can lead to gravitaiional 
memory effects in BSs. JBD theory is closely related to low energy limits of superstring theories [64] which imply 
the primordial production of scalar fields. Then, relics like the dilaton, the axion, combined as axidilaton, or other 
moduli fields could remain in our present epoch as candidates. 

Already the first two papers on BSs provided the two main directions: The history of these hypothetical stars starts 
in 1968 with the work of Kaup [159] using a complex massive scalar field with gravitational interactions in a semi- 
classical manner. The energy-momentum tensor is calculated classically providing the source for gravitation; a very 
detailed investigation of the solution classes was done in [98]. However, one year later, Ruffini and Bonazzola [237] used 
field quantisation of a real scalar field and considered the ground state of N particles. The vacuum expectation value 
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of the field operators yield the same energy-momentum tensor and thus, not surprisingly, the same field equations. 
The two different physical constituents, complex versus quantised real scalar field BS, yield the same macroscopical 
results. It should be noted, however, that the gravitational field g^i, is kept classical due to the non-renormalisability 
of standard GR, or, alternatively, treated it as principal low energy part of some renormalisable superstring model. 

More recently, a BS using a field quantised complex scalar field has been constructed. BSs consisting of pure semi- 
classical real scalars cannot exist because their static solutions in flat spacetime are unstable due to Derrick's theorem 
[80], solutions may arise only if the real scalar field possesses a time-dependence leading to non-static oscillating BSs. 
However, if a fermion star is present as well, a real scalar component can be added and if it interacts with the fermions, 
then a combined boson- fermion star is the result as in the first calculation by T.D. Lee and Pang [182]; cf. Section 
VA. 

A challenge for the BS model is that, so far, no fundamental scalar particle has been detected with certainty in the 
laboratory. Several are proposed by theory: The Higgs particle h is a necessary ingredient of the standard model. 
However, the possible discovery of the Higgs boson of mass Wh = 114.5 GeV/c^ at the Large Electron Positron (LEP) 
collider at CERN [3,20] gives the BS strand of investigation a fresh impetus. 

For the BS, we need a scalar particle which does not decay; or if it decays (like the Higgs, e.g.), one has to assume 
that, in the gravitational binding, the inverse process is efficient enough for an equilibrium, as is the case for the 
/3-decay inside a neutron star. In the latter case, the direct physical predecessor of that kind of a BS is not clear; 
but, of course, unstable Higgs particles could not have formed a BS by themselves. In order to explore that unknown 
particle regime, one can play with the model parameters such as particle mass or interaction constants. Then, BSs 
of rather different sizes can occur: it could be just a 'gravitational atom'; it could be as massive as the presumed BH 
in the central part of a galaxy; or it could be an alternative explanation for parts of the dark matter in the halo of 
galaxies. 

In 1995, experiments [7] proved the existence of the fifth possible state of matter, the Bose-Einstein condensate 
(BEG); in 2001, the Nobel prize was awarded for its experimental realization in traps. BSs, if they exist, would be 
an astrophysical realization with a self-generated gravitational confinement, cf. [158,21]. 

Let us sum up some of the properties of complex scalar field BSs (following an earlier version of [296]) in Table I 
as we shall discuss in Section HI. 

TABLE I. Overview of some complex scalar field BS properties. 



Property 


BS 


Constituents 


Scalars (Bose-Einstein-Condensation) 


Pressure support 


Heisenberg's uncertainty relation 


Size 


Gigantic up to very compact (few Schwarzschild radii): Table IV, Fig. 3 


Surface 


Atmosphere 


Appearance 


Transparent (if only gravitationally interacting) 


Structure 


Einstein-Klein-Gordon equation 


Gravitational potential 


Newtonian weak up to highly relativistic 


Last stable orbit 


None 


Rotation 


Difl[erentially; discrete 


Avoiding a baryonic BH by 


Jet, particle dynamics 


Avoiding a scalar BH by 


Stability, evolution 


Gravitational redshift 


Comparable to neutron star values, but larger due to transparency 


Gravitational (micro-)lensing 


Extreme deflection angles possible (Fig. 2); MACHOS 


Luminosity 


Larger than luminosity of a BH 


Star disruption (tidal radius) 


Yes 


Distinctive observational 


Broadening of emission lines 


signatures 


Gravitational waves 




Cerenkov radiation 
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A. Fundamental scalar fields in Nature? 



The physical nature of the spin-O-particles out of which the BS is presumed to consist, is still an open issue. Until 

now, no fundamental elementary scalar particle has been detected with certainty in accelerator experiments, which 
could serve as the main constituent of the BS. In the theory of Glashow, Weinberg, and Salam, a Higgs boson doublet 
($+, and its anti-doublet ($~, $'^) are necessary ingredients in order to generate masses for the 'heavy photons', 
i.e. the and gauge vector bosons [233]. After symmetry breaking, only one real scalar particle, the Higgs 
particle h := (4>° + $°)/V2, remains free and occurs in the state of a constant scalar field background [229]. As it 
is indicated by the rather heavy top quark [1] of 176 GeV/c^, the mass of the Higgs particle is expected to be below 
1000 GeV/c^. This is supported by the possible discovery of the Higgs boson of mass Wh = 114.5 GeV/c^ at the 
Large Electron Positron (LEP) collider at CERN [3,20]. In order to stabilise such a light Higgs against quantum 
fluctuations, a supersymmetric extension of the standard model is desirable. However, then it will be accompanied by 
additional heavy Higgs fields H^, A^, and a charged doublet in the mass range of 100 GeV/c^ to 1 TeV/c^. For 
the unlike case of a Higgs mass rrih above 1.2 TeV/c^, the self-interaction J7($) of the Higgs field is so large that any 
perturbative approach of the standard model becomes unreliable. Tentatively, a conformal extension of the standard 
model with gravity included has been analysed, cf. [227,126]. Future high-energy experiments at the LHC at Cern 
should reveal if Higgs particles really exist in Nature. 

As free particles, the Higgs boson of Glashow- Weinberg-Salam theory is unstable, e.g. with respect to the decays 
h -I- W~ and h ^ + Z^, ifit is heavier than the gauge bosons. In a compact object like the BS, these decay 

channels are expected to be in partial equilibrium with the inverse processes Z^ + Z'^ ^ h or, via virtual W triangle 
graph, Z'^ ^ h + j, for instance [50], by utilising gravitational binding energy. This is presumably in full analogy with 
the neutron star [136,273,306] or quark star [160,110], where one finds an equilibrium of f3- and inverse /3-decay or 
of a quark-gluon plasma (with similar decay and inverse-decay processes) and thus stability of the macroscopic star 
with respect to radioactive decay. 

The known spin-O-particles in particle physics, although effectively described by a Klein-Gordon (KG) equation, 
are not elementary, they consist of quarks having spin 1/2. Nevertheless, they possess some physical properties which 
can help us to understand which kind of particles can occur in the different BS types. 

The electrically charged pions tt^ with a mass of about 140 MeV/c^ are a typical example for complex scalar 
fields [138,194]. The electrically neutral K°, K° mesons are described by a complex KG field as well; K° and K" are 
distinguished by their hypercharge Y := B + S = ±1, where B = is the baryon number and S is the strangeness 
corresponding to a global U{1) symmetry. In both cases, particles (e.g. 7r+, K*^) and anti-particles (e.g. n~ , K^) occur. 
Generally, Noether's theorem leads to a conserved charge 

Q = e{N+-N-) (1.1) 

related to the number N of particles and anti-particles, respectively, where e defines the absolute value of the charge. 

Let us stress that if the KG equation for a complex scalar field admits a global U{1) symmetry, as in the case for the 
K*\ K'^ mesons, "merely" a conserved particle number N arises and the charge is the total strangeness S. Regarding 
the high degree of instability of pions, K mesons, and Higgs fields, the effect of gravitational binding is essential to 
see if these particles had enough time to form a BS. 

In string theories, there may exist several fundamental scalar particles having different global or local charges. For 
example, for complex scalar particles with global U{1) symmetry or for real scalar particles, the BS charge shall be 
called just boson number N. For complex scalar particles with local U{1) symmetry, the BS has a charge given by 
(1.1). In a BS with temperature near zero, we will assume that all constituents are particles, i.e. no anti-particles 
present, or vice versa. 

In his model, Kaup [159] used a complex scalar field with global U{1) symmetry so that no gauge bosons are 
present. Ruffini and Bonazzola [237], however, took a real scalar field for the general relativistic BS and a complex 
scalar field with global U{1) symmetry for the Newtonian BS description. In the case that U{1) symmetry is local, 
the complex scalar field couples to the gauge field contribution, cf. [138], p. 31; Jetzer and van der Bij [147] called 
their hypothetical astrophysical object charged boson star (CBS). Within the Glashow- Weinberg-Salam theory before 
symmetry-breaking of the group SU{2) x U{1), this U{1) charge describes the weak hypercharge, and not an electric 
one. During the earliest stages of the universe, a complex scalar field with different kinds of U{1) charge (than electric 
or hyperweak) could have been generated. 

There is also the particle/field classification with respect to the spatial refiection P. Both a real and a complex 
scalar field can be either a scalar or a pseudo-scalar. 

Thus, particle characterisation will lead to different consequences for the detection of each BS model. 
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B. Boson star as a self-gravitating Bose-Einstein condensate 

Since Einstein and Bose it is well-known that scalar fields represent identical particles which can occupy the same 

ground state. Such a Bose-Einstein condensate (BEC) has been experimentally realized in 1995 for cold atoms of 
even number of electrons, protons, and neutrons, see Anglin and Ketterle [7] for a recent review. In the mean-field 
ansatz, the interaction of the atoms in a dilute gas is approximated by the effective potential 

umu = jm'. (1.2) 

This leads to a nonlinear Schrodinger equation for ^, in this context known as Gross-Pitaevskii equation. In a 
microscopic approach, one introduces bosonic creation and annihilation operators and b, respectively, satisfying 

[6,6t] = i (1.3) 

and finds that every number conserving normal ordered correlation function (&i • • • &„) splits into the sum of all 
possible products of contractions {bibj) as in Wick's theorem of quantum field theory (QFT). For n = 1 one recovers 
the Gross-Pitaevskii equation, whereas the next order leads to the Hartree-Fock-Bogoliubov equations, see [165] 
for details. Therefore, it is gratifying to note that BSs with repulsive self-interaction C/(|$|) considered already in 
Refs. [210,58] have their counterparts in the effective potential (1.2) of BEC. Thus, some authors [158,21] advocate 
consideration of a cold BS as a self-gravitating BEC on an astrophysical scale. 

Recently, the so-called vortices, collective excitations of BECs with angular momentum in the direction of the 
vortex axis z, have been predicted [313] and then experimentally prepared [192]. Since, quantum-mechanically, 
the z-component Jz = aNh of the total angular momentum is necessarily quantised by the azimuthal quantum 
number a, it is not possible to "continuously" deform this state to the ground state, and by this circumstance, 
contributing to its (meta-) stability. In 1996, axisymmetric solutions of the Einstein-KG equations have been found 
by us [253,211,257,260] which are rotating BSs (differentially due to the frame-dragging in curved spacetime) and 
exhibit, as a collective state, the same relation J = aN for the total angular momentum as in the case of the vortices 
of an BEC and a meta-stability against the decay into the ground state as well; cf. Section II I. 



II. COMPLEX SCALAR FIELD BOSON STARS 



In this Section, in a nutshell, BSs are localised solutions of the coupled system of Einstein and general relativistic 
Klein-Gordon (KG) equations of a complex scalar Depending on the self-interaction potential J7(|$p), BSs have 
received different labels: In 1968, Kaup [159] referred to the localised solution of the linear KG equation as a Klein- 
Gordon geon, whereas in 1987, the same compact object was called mini-soliton star [98] within a series of publications 
[180,181,98,99] where also the terms scalar soliton star or just soliton star were introduced which nowadays are coined 
non-topological soliton star. Before, in 1986, Colpi et al. [58] baptised a self-gravitating scalar field with ]<I>]^ self- 
interaction a boson star. So, in 1989, in a paper on stability, the mini-soliton star was renamed to mini-boson star 
[183]. An additional U{\) charge led to the label charged BS [147]. In [47], different BS models were summarised 
under the label scalar stars. Let us stress that already Ruffini and Bonazzola [237] envisioned a star by calling their 
investigation systems of self- gravitating particles in GR. 

In his perspective paper, Kaup [159] had studied for the first time the full general relativistic coupling of a linear 
complex KG field to gravity for a localised configuration. It is already realized that neither a Schwarzschild type 
event horizon nor an initial singularity occurs in the corresponding numerical solutions. Moreover, the problem of the 
stability of the resulting BS with respect to radial perturbations is treated. It is shown that such objects are resistant 
to radial gravitational collapse for a total mass M below some critical mass and low central density (related works 
include Refs. [73,91,285]). The resulting configuration is a macroscopic coherent state for which the KG field can be 
treated as a semi-classical field. 

The Lagrangian density of gravitationally coupled complex scalar field $ reads 

jCbs = {r + «[9''^(a^$*)(a.$) - Um^)] } . (2.1) 

Here n = 8itG is the gravitational constant in natural units, g the determinant of the metric gfj,^, R ■= g^^ R^v = 
g^'^ i^d^T — d^T^y'^ + T^fj"-Tav'^ — F^i^^FaCT'^) the curvature scalar with Tolman's sign convention [290]. Using the 
principle of variation, one finds the coupled Einstein-Klein-Gordon equations 
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□ ■ 



dU 



□ 



d I $ |2 

dU 

d I $ |2 



$ = 
$* = 



where 



is the stress-energy tensor and □ 
The stationarity ansatz 



(l/VI 9 \ j (^fj. (Vl 9 Ig'^'^dSj the generahy covariant d'Alembertian. 



(2.2) 
(2.3) 
(2.4) 

(2.5) 



(2.6) 



and its complex conjugate describe a spherically symmetric bound state of scalar fields with positive or negative 
frequency uj, respectively. It ensures that the BS spacetime remains static. (The case of a real scalar field can readily 
be accommodated in this formalism, but requires u) = due to $ = $*.) Equations (2.3) and (2.4) give the same 
differential equation. BSs composed of complex scalars are symmetric under particle anti-particle conjugation C, i.e., 
there exists a degeneracy under the sign change lu —to, N —N, e — e, cf. Eq. (2.18), where e denotes the U{1) 
charge and N := J j^d^x the particle number; cf. Section II C. Hence, a BS in total consists of either particles or 
anti-particles but the classical description in field equations above give us an indication of its constituents by the sign 
of N . This persists also in field quantisation where (1.1) holds; cf. Sections IIR and IVA. 
In the case of spherical symmetry, the general line-element can be written down as follows 



ds' 



=A(r) 



-I- I 



d6'2 -I- sin2 



9dip^^ 



(2.7) 



in which the metric is static and the functions ly = iy{r) and A = A(r) depend only on the Schwarzschild type radial 
coordinate r. Prom (2.2)-(2.3), a system of three coupled nonlinear equations for the radial parts of the scalar and 

the (strong) gravitational tensor field results which has to be solved numerically. 

Similarly as in the case of a prescribed Schwarzschild background [78] , the self-generated spacetime curvature affects 
from the KG equation (2.3) the resulting radial Schrddinger equation 



[a„2 - Feff (r*) + a;2 - m2] P = 



(2.8) 



for the radial function P(r) := $e*'^* essentially via an effective gravitational potential Vctf{r*) = dU / {d\^\'^) + 
e^lil + 1) /r2 + {v' - X')e"-^/2r, when written in terms of the tortoise coordinate r* := J eS^-"'^/'^ dr, cf. [210] . Then, 
it can be easily realized that localised solutions decrease asymptotically as P{r) ^ (1/r) exp (^—y/m^ — uP' r) in a 
Schwarzschild- type asymptotic background. As first shown by Kaup [159], cf. [237] for the real scalar field case, 
metric and curvature associated with a BS remain completely regular. 



A. Anisotropy of scalar matter 

That the stress-energy tensor (2.5) of a BS, unlike a classical fluid, is in general anisotropic was already noticed by 
Kaup [159]. For a spherically symmetric configuration, it becomes diagonal, i.e. Tju''($) = diag {p,—Pr, —p±,—p±) 
with 

p= l(a;2p2e- + P'2e-^ + C/) , 

Pr=p-U, 

P±=Pr- P'^e-^ . (2.9) 

In contrast to a neutron star [122,321,273,277], where the ideal fluid approximation demands the isotropy of the 
pressure, for spherically symmetric BSs there are different stresses Pr and p±_ in radial or tangential directions, 
respectively. Ruffini and Bonazzola [237] display in their Fig. 3 the anisotropy of the different stresses for a quantised 
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real scalar field BS. Since it satisfies the same differential equations, we can use this figure here as well. Gleiser [107] 
introduced the notion oi fractional anisotropy af := {pr —p±)/Pr = P''^e~^ / {p — U) which depends essentially on the 
self-interaction. Furthermore, the contracted Bianchi identity W^Tf/ = is equivalent to the equation 

^Pr = -^' {p+Pr-'^iPr-pl.)^ (2.10) 

of 'hydrostatic' equilibrium for an anisotropic fluid, a generalisation of the Tolman-Oppenheimer-Volkoff equation, see 
Ref. [211]. 

By defining the number density n := uje~^^'^P'^ and the differential pressure := Pr — P^-= P'^e~^, Kaup found 

PA = PAip,Pr,n) = p+Pr - Tl^m^ / {p - Pr) (2-11) 

as an algebraic equation of state, valid only for the linear KG case. However, the exact form of this algebraic relation 
depends on the differential equations (2.2)-(2.3). Consequently, these four thermodynamical variables for the BS show 
a nonlocal behaviour upon perturbations depending on the boundary conditions [159]. Furthermore, Kaup derived 
that radial perturbations have to be nonadiabatic. 

For a spherically symmetric BS in its ground state two different layers of the scalar matter are separated by 
P±{Rc) ~ 0, i.e. a zero of the tangential pressure. Near the centre, p± is positive and, after passing through zero at 
the core radius Rc, it stays negative until radial infinity [247]. The core radius i?c is still inside the BS and contains 
most of the scalar matter. Hence, all three stresses are positive inside the BS core; the boundary layer contains a 
matter distribution with > and p± < 0. 

The choice of the potential U implies a difference in the radial pressure Pr — p — U , and, therefore, may affect the 
physical properties of the BS. In Table II, we shall indicate the main proposals of current literature. 

TABLE II. Complex scalar field BS models distinguished by the scalar self-interaction and 
year of the first publication. 



Compact 
Object 


Self-Interaction C/(|$|2) 
[a, /3, A, A(2„+2), ^0 are constants] 


Year, Publication 


Mini-BS: 
Newtonian BS: 
Self-interacting BS: 
BS: 

Non-topol. Soliton Star: 
General BS: 
Sine-Gordon BS: 
Cosh-Gordon BS: 
Liouville BS: 


C/n = m2 $|2 

C/HKG = m2|$|2-a|$|4 + /3|$|6 
;7csw = m2|<D|2 + A|cI>|V2 

;7NTS=m2|$|2(l-|$|2/<l>2)2 

UhKL = C/CSW + . . . + A(2„+2)|$|'"+2 

UsG = am"^ [sin(7r/2[/3v/|W - 1]) + l] 
UcG = am^ [cosh(/3v^|$|2) - l] 
i/L = am2 [exp(/32|$|2)_i] 


1968, Kaup [159] 

1969, Ruffini-Bonazzola [237] 
1981, Mielke-Scherzer [210] 

1986, Colpi-Shapiro-Wasscrman [58] 

1987, Friedberg-Lee-Pang [99,184] 

1999, Ho-Kim-Lee [131] 

2000, Schunck-Torres [262] 
2000, Schunck-Torres [262] 
2000, Schunck-Torres [262] 



The index of ?7hkg in Table II stands for Heisenberg-Klein-Gordon. As indicated, the scalar potential is a function of 
1$]^. This means that the field $ can be either a scalar or a pseudo-scalar, depending whether for spatial refiections 
there exist an unitary operator V for which V^{x,t)V^ = ±$(— holds. For instance, the axion a is a pseudo- 
scalar. The invariance of the Lagrangian under V depends on whether or not the self-interaction is even or odd. It 
should be noticed that the inclusion of |<I>|^ or higher order terms into the potential implies that the scalar part of 
the theory is no longer rcnormalisable. Fig. 1 shows the Feynman diagram for a l*!']^ interaction. 

In the Mielke-Scherzcr papca- [210], actually the BS is composed from several complex scalars which are in the same 

ground state of a t'Hooft type monopole configuration $^ ~ i?(r) Pj'^' (cos ^), where / indicates the SO{2N) group 
index. In their calculation, an averaged energy- momentum tensor is used as proposed by Power and Wheeler [231], 
p. 488, deriving an angular momentum term in the field equations. 

Regarding the different BS labels, a historical remark is in order. In 1984, Takasugi and Yoshimura [282] calculated 
the gravitational collapse of a cold Bose gas by using the Tolman-Oppcnheimer-Volkoff equation and called their result 
a cold Bose star, in the same year, the same name was given to Ruffini and Bonazzola's BS [40]. The label boson star 
was first introduced by Colpi et al. in 1986 [58]. The reason for that may have been the similarity of the mass units 
with neutron star mass imits which in both cases are the order of magnitude A^fpj/m^ (for details see Section II C). 
For if the scalar field mass m is in the order of GeV/c^, we find for a BS the same mass range as that of a neutron 
star. Due to this argument, Kaup's BS with mass unit Mp^/m seems to be a small object, and it has been called 
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mini-BS [183]. Wc would like to stress that the size and total mass of a BS depends strongly on the unknown value 
of m and on the issue whether there are scalar self-interactions. For corresponding m's, even a so-called mini-BS may 
have masses observed in active galactic nuclei; cf. Sections II M, IIIC, HID, IIIG. The existence of mini-BSs has also 
been proven mathematically [31]. 



B. Charged boson star 



If the complex scalar field is minimally coupled to local U{1) gauge fields, the Lagrangian density necessarily 
(cf. [138]) adopts the form 



.^CBS 



1 

2k' 



R + 



1 



(2.12) 



In accordance with the gauge principle, the coupling between the scalar field and the [/(l)-valued 1-form A = Afj.dx^, 
is introduced via the gauge and, for scalar fields, general covariant derivative -D^$ — 9^$ -f- ieA^^, where e denotes 
the U{1) coupling constant. Furthermore, the Maxwell type term for the two-form F := ^Fij,i,dx^ Adx" occurs, where 



dfjiAi, — d^Afj, is the Faraday field strength. 



The expanded form 

(£)^$)* (£>'*$) = (a^$)(a''$*) + ieA''$a^$* - ieA'*$*5^$ + e^A^yl^l^l^ 



(2.13) 



of the kinetic term in (2.12) helps us to understand the gauge interactions in flat space-time for weak couplings. In 
Fig. 1, the two Feynman diagrams in order e and order in the U(l) gauge coupling are illustrated for flat spacetime. 






FIG. 1. Feynman graph for the A|<1>|* sclf-intcraction in Minkowski spacetime [left]. Feynman graphs for the eA*^^* (or 
eA^$a^$*, resp.) [middle] and for the e^^^A^|$|^ [right] self-interaction. 



Semi-classical solutions for (2.12) have been found by Jetzer and van dcr Bij [147] for an electrically charged 
complex scalar field. These spherically symmetric configurations, called charged BSs (CBSs) [147], can only exist if 
the gravitational attraction is larger than the Coulomb repulsion, i.e. if a = e^/47r < e^^j^ = m^/Mp^, where a is the 
fine structure constant; cf. the stability investigation in [149]. A macroscopic CBS would accrete matter of opposite 
charge very rapidly or be destroyed if the accreted matter consists of anti-particles [185]. There may be vacuum 
instabilities induced by a CBS [154] as we discuss now. 

For particle-like CBSs, a screening mechanism can occur [154]. If the CBS is more compact than the Compton 
wavelength of the electron and supercritically charged, then pair production of (i) electrons and positrons or of (ii) 
pions occurs. For fermionic pair-production, the critical CBS charge is about Zc ^ 1/a, whereas for bosonic pairs, 
it is found that Zc ~ 0.5/a. If a CBS has a supercritical charge, the particles with opposing charge from the pair- 
production screen every surplus above the critical charge whereas their partners leave the newly formed compound 
object. In [154] an example of a CBS is given with 66 particles very concentrated at the centre, and a cloud of 
10 screening particles, e.g. pions, almost outside of the CBS (m^ — lOOOw^); the screened CBS has a net charge 
of 56e. There is also almost no influence on the physical properties of CBS from the condensate (apart from its 
electromagnetic properties at large distances). For CBSs with constituent mass m comparable to TOtt, the vacuum is 
no longer overcritical, and no pair production occurs. Furthermore, the pair production for the bosonic sector is three 
order of magnitudes faster than for the fermionic one, so that a pion condensate will most likely be created first. The 
CBS radius in this example is less than 0.1 fermi. Fiirthermore, it is a speculation that CHAMPs (CHArged Massive 
Particles) could actually form CBSs, if such charged particles like pyrgons or maximons would exist [154]. 

A first-order perturbative approach of a charged bosonic cloud has recently been examined [71] by which analytical 
formulae for mass, charge, and radius are achieved if we have oj = m, i.e., the cloud is in an energy state even above 
the ones of excited CBS states. This investigation is the starting point of the dynamics of CBS formation as continued 
in Ref. [72]; cf. Section IIIF. 

For the detection of CBSs, we should distinguish the different meanings of a U{1) symmetry. For example, within 
the standard Glashow-Weinberg-Salam theory, before symmetry-breaking of the group SU{2) x U{1), the U{1) group 
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describes the weak hypercharge. Hence, during the earliest stages of the universe, a complex scalar field with different 
kinds of U{1) charge (than electric or hyperweak) could have been generated. If we allow for a spontaneous symmetry 
breaking, the Higgs field in the standard model becomes a real neutral scalar field afterwards, and possible mixing 
terms are eliminated in the unitary gauge. Since, for a bound state like a BS, the weak coupling expansion does not 
apply, in principle lattice field theory comes into force, cf. [50]. 

In a pre-CBS topic, stationary axisymmetric solutions are given in [19]. In the physically different but mathe- 
matically similar context of the Ginzburg-Landau system, cylindrically symmetric solutions for charged scalar fields 
without self-gravitation were found in [221,222]. 

C. Critical masses of boson stars 

There is an essential distinction between real and complex scalar fields in the linear KG equation: 

• The Lagrangian of a real scalar field is invariant only under the discrete symmetry $ ^ $' = — $, 

• whereas the dynamics of a complex scalar $ := ($i + i^2)/V^ is invariant under the global 0(2) ~ U{1) 
symmetry $ $e*", where a is a constant. 

The Noether theorem associates with each contimioiis symmetry a locally conserved current S^j^ = and "charge" 
(1.1) which commutes with the Hamiltonian, i.e. [Q, H] = 0. For a complex scalar field, the first "constant of motion" 
of our coupled system of equations is given by the invariance of the Lagrangian density (2.1) under a global U{1) 
transformation $ — > $e'" and so, constraining the form of the self-interaction potential. From the associated Noether 
current j^, the particle number N arises: 

f = ^VW] 9''''{^*d.<^ - '^d,^*] , N:= j fd^x ; (2.14) 

cf. [111,112]. The second "constant of motion" arises from the Abelian group of time translations, which gives rise 
to the concept of energy, cf. [209]. In a rest frame, this is proportional to the total gravitational mass M of localised 
solutions. Since these are asymptotically isolated configurations we can use Tolman's expression [289] : 




(2.15) 



In the case of a spherically symmetric CBS, where the covector of the electromagnetic potential is given by = 
[C(r), 0,0,0], mass and particle number become 

(2.16) 



g(.+A)/2^2 ^ (2.17) 



(2.18) 

Equivalently, the BS mass satisfies M = Aw pr^dr which enters as a mass parameter in the metric potentials 
gi/(r) g-A(r) ^ X _ 2M/r of the vacuum Schwarzschild solution, or in the asymptotic regions where the BS matter 
tends to zero. 

Since BSs are macroscopic quantum states, they arc prevented from complete gravitational collapse; by the Heisen- 
berg uncertainty principle. This providers us also with a crude mass estimate: For a boson to be confined within the 
star of radius i?o, the Compton wavelength of the scalar particle has to satisfy A* = (27r?i/mc) < 2Rq. On the other 
hand, the star's radius should not be much less than the last stable Kepler orbit 3i?s around a BH of Schwarzschild 
radius Rs ■= 2GM/c^ in order to avoid an instability against complete gravitational collapse. 

For the mmi-BS, let us therefore assume an effective radius of Rq = (7r/2)^J?s = 2.47i?s- So far, the factor (7r/2)^ 
has no theoretical foundation, but then the estimate 



oo 

M = '^'K j {p + pr + 2p±)e'-''+^'>/^r^dr 



oo 

= 4./ 



OO 

N = 4n j P^{iJ + eC)e(^-'^)/V2 dr . 
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Merit = (2/7r)M|i/m = 0.6366 M^Jm > Mxaup = 0.633 M^Jm , (2.19) 

cf. Ref. [152], provides a rather good upper bound on the Kaup limit MKaup which is the numerically determined 
maximal mass of a stable BS. Here Mpi := y^hc/G is the Planck mass and m the mass of a bosonic particle. As is 
typical for solitonic solutions, the star becomes heavier for weaker coupling C/(|$|), i.e. for smaller constituent mass 
m in a mini-BS. 

In building star-sized BSs, one needs particles of ultra-low mass m or a non-linear [210,58,213] repulsive self- 
interaction ?7(|<i>|)csw = m^|<i>|2 ^1 + A|<i>|2/8^, where A = iX/m'^. In the latter case, for very smah A = 

AMpi/(47rm^), the scalar field scales still as $ ~ Mpi/-\/T67r and it turns out [213] that the Kaup hmit (2.19) 
can be applied again, but for a rescaled mass m — > rriresc := m/ -\/l + A/8. Then, in this regime, the maximal mass of 
a stable BS is approximately 

Merit ^ - + A/8 — . (2.20) 
TT m 

For large A, the scalar field scales as $ ~ Mpi/\/A, so that U (|$|)csw — (1 + 27r)m^Mpj/A; the rescaled mass is now 
m mresc := m^/(T+2^T)/A and we find 



M„H 2 = 0.236^/A ^ ; (2.21) 

TT V 1 + 27r TO TO 

cf. Fig. 2 from Colpi et al. [58] where, as well, the value for the critical mass of 0.22 is derived in Eq. (18). The 
complex scalar field can be redefined as a dimcnsionlcss field a = ■\/k/2 $ — \/4ttG $ and the values of a for stable 
BSs are in the order of unity. Then, depending on the constituent mass and the coupling constant, the critical mass of 
a BS can reach the Chandrasekhar limit Mch := M|[/m^ 2:^ Mq (if m is about the neutron mass), where Mq denotes 
the mass of the sun, and so can possibly even easily extend the limiting mass of 3.23 Mq for rotating neutron stars 
[195]. 

TABLE III. Critical mass and particle number for complex scalar field BS models. 



Compact 


Critical Mass 


Particle Number 


Object 


Merit 


A^crit 


Fermion Star: 


Mch := M|i/to2 


~ (Mpi/to)^ 


Mini-BS: 


MKaup = 0.633 M^i/to 


0.653 (Mpi/to)2 


BS with UcG (/3 = 1): 


M = 0.638 M|i/to 


0.658 (Mpi/to)2 


BS with UsG (/3 = 1): 


M = 0.620 M2j/to 


0.639 {Mpi/m)^ 


Rotating mini-BS: (a = 1) 


1.31M2j/to 


1.38 (Mpi/to)2 


Rotating mini-BS: (a = 2) 


> 2.21 M|i/to 


> 2.40 (Mpi/to)2 


BS: (A > 47rmVM2j 


(l/7r\/2^)\/AM3j/m2 


~ \/A(Mpi/to)3 


Non-topological Soliton Star: 


0.0198(Mi*j/m4>g) 


0.0055(M|i/to2$3) 


CBS: (A = 0) 


~ 0.44/Veerit - e M^^/m 


- 0.44/Veerit - e {M^./mf 


CBS: (A > 47rmVM2j 




~ (Mpi/m)3/2 


~ 0.226/ Veerit - e x/AM^^/m^ 


~ \/A/Vecrit - e (Mpi/m)^ 






~ (Mpi/to)5/2 



In the case of a CBS (Table III), the charg measured in units of order of magnitude Mp\/m [147], 

so that CBSs have actually a fractional mass law [251]. The non-topological soliton stars can easily have very large 
masses; but compare the discussion of fermion soliton stars in [62,63] where the order of magnitude of a neutron star 
can be found, cf. Section V A. 

In general (cf. Table III), a BS is about the huge factor {Mpi/m,) more massive than the mini-BS, as long as 
(Mpi > to). However, this depends also on the value of the coupling constant A, which even for the real Higgs scalar 
of the standard model is experimentally quite unconstrained. In astrophysical terms, the critical mass of a stable BS 
is Merit — 0.127\/AMch = 0.1\/A (GeV/mc^)^ Mq [58]. For cosmologically relevant (invisible) axions (described by a 
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real pseudo-scalar a) of rua ~ 10 ^ eV/c^, an axion star would have the ridiculously large mass of Merit ^ 10"^^ VXMq 
[58]. 

In Table IV, wc show several possible sizes which such localised configurations could have. We consider two models: 
on the left, the BS case with an effectively large self-interaction of [58] where A ^ A-kw?' /M^^, on the right hand side, 
the self-interaction constant A does not fulfil this condition or we have a mini-BS situation, for both cases, we can 
use the values of M and R for A = 1. The order of magnitude of the critical mass and radius are given, hence, R is 
calculated by the Schwarzschild radius IGMjc? and so the BS configuration is effectively general relativistic. Table 
IV provides examples of rather different physical situations. For example, in order to obtain a BS with about a solar 
mass Mq = 10-''° kg = 10'" GeV/c^, a scalar field mass of m - 10"!" eV/c^ for A = 0, or m cx 1 GeV/c^ with A = 0(1) 
is needed; for the latter situation, we find the Colpi et al. condition A » 10"'^'', hence, the self-coupling can only be 
neglected if it is extraordinarily tiny. In this example, the self-interacting scalar particle has a mass comparable to a 
neutron, leading to a BS with the dimensions of a neutron star, cf. Section HID. Furthermore, we recognise that the 
free scalar particle leads to a solar mass BS as well so that it is not so "mini" as its name mini-BS implies. We give 
also values for a mean density p ~ MjR? . If we reduce the self-interacting scalar mass further, to m ^ 1 MeV/c^, 
then we find a localised object that (for A = 1 again) has the size of the Sun but consists of 10^ solar masses; this is 
reminiscent of supermassive black holes, for example as in active galactic nuclei. A mini-BS of the same dimensions 
would need an extremely small scalar mass, which, for the time being, has no observational support. By reducing m 
further, we get a BS which could play a role as a galactic halo, cf. Section II M. 

Let us now increase m. With self-interaction and m = 100 GeV/c^ (about the possible Higgs mass), the BS is just 
about 0.1 m but with a mass of 10^^ kg; for the mini-BS with the same bare mass, we find an ultra tiny atto-mcter 
BS, clearly below one Fermi. On the left hand side of this same line, in order to derive the same object with an 
effectively self-interaction constant, a very heavy scalar particle is needed which is beyond experimental estimates as 
well. Other not mentioned values of m can easily be inter- or extrapolated. 

TABLE IV. Physical sizes of complex scalar field BS. 



effective SI 
Ucsw ^ 


BS configuration 

■*— needs — > 


(takt 


uneffective SI or 

C/k 7^0 
i A = 1 in middle row) 


m 


= 10^"-'' GeV/c2 


M s 


a 10'5^/A kg, R w 


l\f\ atto m. 




i 1063/A kg/m3 


TO = 


100 GeV/c2 


m 


= 106-5 GeV/c2 


M s 


^ 10"VA kg, R R 


i iVa A, 




i 10^7 A kg/m^ 


TO = 


1 keV/c2 


m 


= 10^ GeV/c^ 


Mi= 


i 10^0 \/A kg, R R 


i 100\/A nm, 




i IO^VA kg/m^ 


m = 


1 eV/c^ 


m 


= 100 GeV/c2 


Mf= 


i W-^VXMq, R 


K. 0.1 V\ m, 




i IO^VA kg/m3 


TO 


10-6 eV/c2 


m 


= 1 GeV/c2 


M s 


^ 10° VxMq, r 


10 \/A km. 




i lO^VA kg/m3 


■ = 


10-1° eV/c2 


m 


= 1 MeV/c2 


M s 


a 10^ Vxmq, R', 


« 10^ VA km. 




i 10^ /A kg/m3 


TO = 


10-16 eV/c2 


m 


= 1 keV/c^ 


M s 


s IOI^TaMq, R 


=s IO-^VA pc, 




i 10° /A kg/m3 


TO = 


10-22 eV/c2 



A different approach to estimate the mass of a mini-BS is already described in [152]. There, analytic bounds are 
derived on the ground state of a system of A'' self-gravitating and relativistic bosons in Newtonian approximation, 
without distinguishing between real or complex scalars. The Hamiltonian for an assembly of bosons is 



1=1 i>j ' ' 

as implicitly used in Chandrasekhar's calculation for white dwarfs [273]. The improved estimation in Ref. [234] for 
the lower and upper bounds on the maximal mini-BS mass is given by 

0.8468 M|i/to <M < 1.439 M|i/to . (2.23) 

Since in the Hamiltonian general relativistic effects are ignored, the bounds are above the Kaup limit Mxaup. 

D. Effective radius 

The surface and thus the radius of a neutron star is characterised by a vanishing (radial) pressure [223,273]. In 
contradistinction to neutron stars, due to the exponentially decreasing scalar field in a BS, we can merely speak of 
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an exosphere of a BS which corresponds to the highest layer of a planet; alternatively, we can call it the halo of a BS. 
For an effective radius of a BS, several proposals have been delivered in the literature. 

The total mass M can be used to define a radius where a certain percentage of mass is reached: R = M(95%) or 
even R = M(99.9%), e.g. [65]. Alternatively, the energy-density p is used in 

J^pr^dr 

Rl ■— -TOO 53- ) [2.24) 

Jo Pr^df 

or with help of Tolman's expression (2.15) [98]. Taken the close relation to quantum mechanical wave fimction into 
account, Mielke [204] proposed the square root of the expectation value of the space operator as BS radius 

^2-<N >- J--|$|2rf, ' - |ci>|2d, ' (2.25) 

the latter being even simpler. All three definitions were investigated numerically in [247] for spherically symmetric 
BSs with the result that Ri < R2 < R3 for one BS state; with increasing cent ral d ensities the BS radius decreases, 
there is a minimum at P{Q) ~ 0.95-\/2/«; and a maximum around P(0) 1.5-^/2/k. However, both these values are 
deep in the unstable BS states. 

Applying the conserved current (2.14), another definition 

1 r°° 

Res:=^ fr^dr (2.26) 



was introduced in the case of CBSs, cf. [147]. The reason is the long range nature of the electromagnetic force; in 
this case, the definition i?i would yield an infinite radius. Let us notice that the Schwarzschild metric gives a general 
lower bound on the radius of any type of object, namely R > M/Mp^. 

In all these definitions, the dimensions of the normalisation and the integrated densities cancel each other, except 
for the radial variable r. Therefore, the effective radius of a mini-BS is proportional to the Compton wavelength 
A$ = 2TTh/mc of the bosonic constituent. This applies also to the case of a BS, except that we approximately can 
use the rescaled mass, leading to A$^^^^ ~ \/A(27r ?i/mc), which, for large A, is about (Mpi/m) (a huge factor) times 
larger. Thus, the corresponding effective radius R^s justifies the distinction between a mini-BS and a BS. 

The situation for a non-topological soliton star is different in this way that it contains an interior part where the 
complex scalar field is about the constant false vacuum value ^>o, a shell of width 1/m over which the complex scalar 
field decreases to zero, and the exterior vacuum. The radius is then defined by the position of the shell. 



E. Stability and Heisenberg's uncertainty principle 

In the case of a "normal" star, the equation of state is approximately that of an ideal gas, i.e. the ratio of pressure 
p to density p is a function of temperature. For a fermion star, as exemplified by white dwarfs or neutron stars, p/p 
is independent of the temperature; for such degenerated matter, the high pressure of the electron or neutron gas, 
respectively, originates from Pauli's exclusion principle, resulting in a very high Fermi energy. 

Bosons are not susceptible to the exclusion principle; instead, they can all occupy the same ground state. That 
bosons cannot be localised within their Compton wavelength is guaranteed by Heisenberg's uncertainty principle 
[107,186]. 

In atoms as well as in a BS, the same principle applies and is the basic warranty for its stability against collapse. 
Prom a macroscopic point of view, our mathematical model shows that the macroscopic scalar radial pressure works 
against gravity. In a nutshell, a cold BS is a huge Bose-Einstein condensate, as we demonstrated in more detail in 
Section IB. 



F. Stability, catastrophe theory, and critical phenomena 

According to Derrick's theorem [80,5], there are no stable time- independent solutions of scalar fields in flat 4D 
spacetime. One ingredient in order to circumvent this no-go theorem is to allow for a time- dependent phase factor 
exp(— icjt) in the stationarity ansatz (2.6) for complex scalars. For BSs, another factor is the attraction of the 
gravitational field which provides a counter-balance to the effective repulsion due to kinetic energy and corresponding 
parts in the self-interaction U ($) of the scalar field. 
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Moreover, for such soliton-type configurations kept together by their self-generated gravitational field, the issue of 
stability against gravitational collapse is crucial. In the spherically symmetric case, it was shown by Gleiser [107,109], 
Jetzer [148 151], as well as Lee and Pang [183] that BSs with small central densities and masses below the Kaup 
limit are stable against small radial perturbations. (Kaup [159] investigated the stability with respect to radial 
perturbations as well but came to a different answer.) Shortly afterwards, we have established via catastrophe theory 
[170-174,248] an even simpler way to prove stability of star solutions in general. The problem of non-radial pidsations 
of a BS was mathematically formulated [166] and later quasinormal modes calculated [316]; cf. Section IIIF. 

Numerical simulations [123,124] of a mini-BS together with a surrounding spherical shell consisting of a massless real 
scalar field have also led to qualitative stability predictions, comparable to results in [171] as we show below. In the 
simulation, the real scalar field implodes towards the mini-BS centre; both fields are merely interacting gravitationally. 
The pulse then passes through the origin, explodes and continues to infinity. Meanwhile, the mini-BS is compressed 
into a state which ultimately forms either a BH or disperses, depending on initial data. By varying the initial 
amplitude of the real scalar field, a critical solution is defined as lying between initial data resulting in BH formation 
and initial data giving rise to dispersal; cf. the numerical work on gravitational collapse of a massless real scalar in 
[113]. 

As starting point of the numerical solution, a stable mini-BS solution is taken, and perturbed by the contribution 
of the real scalar field where two features are important: (i) the initial amplitude, the so-called critical parameter p, 
and (ii) the form of initial data (gaussian or kink). Because of the choice of perturbation and a parameter value p 
close to a critical value p*, the stable mini-BS is transformed into an intermediate oscillating unstable mini-BS state, 
the so-called critical solution; cf. [54]. The closer p is to p* , the more time the unstable mini-BS oscillates; but, for 
p < p*, the mini-BS always "explodes" eventually, i.e. disperses all matter to infinity, leaving just a diffuse remnant 
with low mass. For p > p*, the perturbed mini-BS forms a BH afterwards, where a finite minimal BH mass is found. 
Due to this and a scaling law relating the lifetime r of a near-critical solution to the proximity of the solution to the 
critical point [r ~ — \n{p — p*)], the critical solutions belong to a so-called type I class. 

The critical mini-BS solutions are a hybrid of an unstable mini-BS solution and a small halo near the outer edge 
of the BS; excited mini-BS states are not the explanation, they do differ significantly. The halo part seems to be just 
an artifact of the collision with the real scalar field. 

Additionally in [123] , the numerical simulation method is used to investigate and verify earlier results on the stability 
against radial perturbations of mini-BS states with small and large central densities. The numerical simulation starts 
with a stable mini-BS solution. If the critical mini-BS solution is reached, and if one ignores the matter in its halo 
part, the rest is described by an unstable mini-BS solution with less mass than the stable star had, but larger central 
density value (see Fig. 16 in Ref. [123]). In fact, the numerical method of the linear perturbation theory of [109] 
is compared with the simulation method. Both results coincide with each other, but, moreover, explicit values of 
the frequencies of the fundamental mode and the first harmonic mode are given. For every existing mode, there 
are cither two different stable mini-BS states, or just one unstable one. Because the mode frequencies may have 
observational consequences, we mention them in Section IIIF. Here the binding energy M — mN has no meaning 
[123] in understanding the transitions in the BS stability or the dynamics of the simulations. 

These new simulations are related to the results of an earlier paper [171] where the stability was described qual- 
itatively by catastrophe theory. In numerical simulations, instability arises when the square of the mode frequency 
is becoming negative, i.e. there exist an exponentially increasing mode, and, by this, the mini-BS is destroyed. This 
starts to happen at extrema in the diagrams (M, P(0)) or {N, P{0)), cf. Fig. 3 in this review; these extrema belong to 
so-called cusps in the diagram (Af , A^) , also called bifurcation diagram. By both the numerical investigation and the 
catastrophe theory, one finds that, on the first cusp, one stable mode becomes unstable; on the second cusp, a second 
mode becomes unstable, and so on. In catastrophe theory, the appearance of each cusp is correlated to one Whitney 
surface (for one corresponding mode) embedded in a higher dimensional parameter space. The bifurcation diagram 
[M,N) arises from a particular projection of this high dimensional parameter space. In Fig. 2 of [171] which is the 
alternative bifurcation diagram (M — mN, N) using the binding energy instead of the BS mass, three cusps are shown 
which are labeled A^ci < A^C2 < A^cs; A^C3 is the maximal particle number corresponding to the onset of instability; 
A^ci belongs to the second cusp (or the next minimum following Nc3, resp.), and Nc2 to the third cusp (or the next 
maximum following Nci, resp.). In Section VI. C of [171], the physical behaviour of non-equilibrium BSs is considered 
for Nc2 < N < Nc3- There, for constant BS mass, two BS states exist, one stable and one unstable. Perturbed BSs 
will exhibit oscillations around these two BS states before such a state either collapses to a BH or disperses. This 
depends on the strength of the perturbation, hence, on the parameter p approaching its critical value from above. 

In [123], numerical simulations for exactly this particle number regime are carried out, and we see good agreements 
between the simulations and the general predictions from catastrophe theory. It would be interesting to see whether 
the numerical simulations for the other regimes could be reconfirmed: (a) for A^ < A^cij only oscillating solutions (no 
collapse, no dispersal), (b) for A^ci < N < Nc2, a pulsation behaviour of perturbed BS states consisting of at least 
two different oscillation modes. 
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G. Newtonian approximation - Gravitational atom I 



In a nut-shell, a BS is a stationary solution of a (non-linear) KG equation in its own gravitational field; cf. [203,206]. 

We treat this problem in a semi- classical manner, because effects of the quantised gravitational field are neglected. 
Therefore, a Newtonian BS was also called a gravitational atom [92]; hence, the excited states of a BS is the topic 
of this subsection. The results of the following subsection show that a general relativistic BS (i.e. one with strong 
self-gravity) as well can be designated as 'atom'. Thus, our 'gravitational atoms' represent coherent quantum states, 
which nevertheless can have macroscopic size and large masses. The gravitational field is self-generated via the energy- 
momentum tensor, but remains completely classical, whereas the complex scalar fields are treated to some extent as 
Schrodinger wave functions, which in quantum field theory are referred to as semi-classical as well. 

Since a free KG equation for a complex scalar field is a relativistic generalisation of the Schrodinger equation, we 
consider for the ground state a generalisation of the wave function [127,207] 

\N, n, I, a >: $„;„(t, r, 0, ^) = RUr) Y{\e, (^)e-'('-"'/'*)* 

^ i?f(r)Pj"(cos6')e-*"*'e-'('"'"/^^* , (2.27) 



4-77 

where R^{r) gives the radial distribution, Yi°'{9, ip) are the spherical harmonics, PI^{cos9) are the normalised Leg- 
endre polynomials, and \a\ < I are the quantum numbers of azimuthal and angular momentum. Each mode has its 
characteristic frequency Wni- As is well-known, in the Newtonian limit Einstein's equation reduces to the Poisson 
equation 

= -K.p/2 , (2.28) 

where V is Newton's gravitational potential, p = A''m|$p is the matter density, and N the particle number. Under 
the condition that all bosons are in the same state, the KG equation reduces to 

A$„,„ + 2m{uni + mV)^nia = , (2.29) 

hence, a Schrodinger- type equation. The relation for the conserved particle number N = m J \ ^nia\^r'^ sin OdrdOdip 
gives us a normalisation condition on the wave function R"{r). No self-interaction of the scalar field is introduced. 

The coupled system of (2.28) and (2.29) has to be solved numerically. A factorization as in (2.27) of the potential 
V and of the scalar field $nia yields two coupled ordinary differential equations for the radial coordinate. 

Several solutions have been derived so far. The nodeless ground state with (n = 1,Z = a = 0) was calculated first 
by [237] and can also be found in [92] whereas in [98] even Newtonian solutions up to 5 nodes are obtained. In Fig. 3 
of [92], we can see very nicely the difference in the Newtonian and GR calculation. Additionally in [92], a BS was 
calculated where scalar matter is distributed within the ground state plus some matter in the lowest mode with nonzero 
quadrupole moment (the 3d mode with n = 3 and Z = 2) so that gravitational waves can be radiated; cf. Section 
IIIF. In [274], non- vanishing azimuthal angular momentum a is included, and so, we have rotating Newtonian BSs; 
explicitly, solutions for the combinations {I = l,a = 0), {I = l,a= 1), (/ = 2, a = 0) are presented. 

In cylindrical coordinates {r,z), the system (2.28), (2.29) is an elliptic boundary eigenvalue problem, and a finite 
element method is used for solving the problem numerically in [265] . Due to the invariance of the differential equations 
with respect to z ^ —z, V ^ V,^ ^ ±$, there are positive and negative parity classes where for the first the 
spherically symmetric solution is derived. Solutions with negative parity show two equatorial symmetric blobs of 
scalar matter above and below the equatorial plane, whereas the scalar field is equatorially anti-symmetric. This 
is comparable to a hydrogen atom in the (n = 2, ? = 1, a = 0) state, cf. [37]. Whereas this solution describes an 
axisymmetric Newtonian BS, also non-axisymmetric Newtonian BSs were obtained [319]. There, the scalar field is 
anti-symmetric with respect to one of the following parity transformations: (A) {x,y,z) {x,y,—z), {x,—y,z); (B) 
{x,y,z) {x,y,—z), (x, — j/, z), {—x,y,z). The type A solutions show a BS with four blobs with maximal scalar 
matter density in the {x = 0)-plane. The type B solutions possess eight blobs of scalar matter with zero density in 
the {x = 0)-plane; in both cases, the density vanishes also for y = and 2 = 0. 

In [317], the general relativistic extension of the axisymmetric mini-BS solutions of [265] were exhibited. These 
negative parity solutions are classified by the derivative $r(^ = 0,^ = 0). The critical mass is given by Merit = 
1.05M2j/m. 

An analytical solution for a BS in the limit A 1 is determined in [178] for the rescaled scalar field 0-*(a;*) = 
^J (70 sin {\/2xst) I \/2"a;, , where 70 is the absolute central value of the Newtonian potential. 
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H. Boson stars - Gravitational atom II 



In the case of the general relativistic BS, several excited BS states have been computed. The first solutions with 
nodes, i.e. principal quantum number n > 1, were determined by Miclkc and Schcrzcr [210] who, motivated by 
Heisenberg's non- linear spinor equation [127,207], added self- interacting terms describing the interaction between the 
bosonic particles in the BS configuration. 

The work of Friedberg et al. [98] investigated solutions with up to 56 nodes within the mini-BS scalar field. Solutions 
with nodes possess a similar behaviour in the mass-particle-number function as solutions without nodes, i.e. a maximal 
mass or cusps in the mass-particle-number diagram were found. For the first cusp, the maximal mass depends linearly 
[98] on the particle number. There as well, the solutions with zero nodes are compared with solutions with ten nodes. 

Recently, phase-shifted BS solutions have been constructed [125]. Normally, the BS complex scalar field $(r, = 
P{r) cos(a;t) + iP{r) cos{Lut + S) has a phase of S = tt/2 between its real and its imaginary part. For S ^ ±tt/2, such 
solutions are investigated and show a BS like behaviour. Long-term numerical evolutions indicate that phase-shifted 
BSs are stable but periodic; mass-energy is exchanged between the real and imaginary fields. It seems that BSs can 
be more generic than previously thought. For 6 = 0, the oscillating BSs are approximated [268]. 

I. Rotating boson stars - Gravitational atom III 

In Refs. [257,261,251,253,211], we proved numerically that regular stationary and axisymmetric solutions, which 
can be regarded as rapidly rotating BSs, exist in GR. In the Einstein-KG system (2.2), (2.3), we use an isotropic 
stationary axisymmetric line element 

ds2 ^ ^^^^2 _ e)dtd<fi - l{r, e)d(p^ - e^^'''^^ {dr^ + r^dO^) , (2.30) 

where f,k,l,fi arc some functions to be determined numerically. In order to find rotating BSs in ground state, the 
scalar field ansatz is 

$(r, 9, if, t) := P(r, 6')e-'"'^e-*"* , (2.31) 

where the uniqueness of the scalar field imdcr a complete rotation <!•(</?) = + 2tt) requires a = 0, ±1, ±2, . . ., as 
in quantum mechanics. For the resulting nonlinear elliptical system, we applied a finite difference method and found 
solutions for the mini-BS with a up to 500. Sequences of rotating mini-BSs with a = 1 and a = 2 were constructed 
in [318] where every BS is characterised by the maximal value of the scalar field which is usually not on the rotation 
axis. For o = 1, the energy density has a non-vanishing value on the rotation axis, whereas, for a > 2, it has to 
vanish identically there, i.e., vacuum predominates; cf. the situation of the rotating BS in (2-|-l)-dimensions in Section 
no. The maximum mass increases above the Kaup limit: for o = 1, it is 1.31Mpj/m and, for a = 2, it is at least 
2.21M|i/m. In [199], rotating BSs were calculated with a numerically faster multigrid method. 
The total angular momentum J of a BS is given by 

J = aN , (2.32) 

i.e. a BS can change its rotational state only in discrete steps resembling gravitational atoms. Hence, the angular 
momentum J is quantised by the azimutlial quantum number a. The lifting of the angular momentum degeneracy, 
familiar from quantum mechanics, is a gravito-magnetic effect [198] due to the rotating frame. It is, however, surprising 
that this still holds for macroscopic configurations. Eventually, it is now understandable why no slowly rotating BS 
states near the spherically symmetric ones could be found [162]. 

If we imagine what happens if a BS starts to rotate, we find out that the BS rearranges its structure from the sphere 
to a torus. The BS in this situation looks like the (n = 2, Z = 1, a = 1) state of a hydrogen atom, cf. [37]. Qualitatively, 
due to the finite velocity of light and the infinite range of the scalar matter within the BS, it follows that the form of 
rotation has to be differential, not uniform. The rotating solutions of the coupled Einstein-KG equations represent the 
field-theoretical pendant of rotating neutron stars which have been studied numerically for various equations of state, 
different approximation schemes [100,60,87] as well as with quadrupole moments [195] as a model for (millisecond) 
pulsars. 

A general insight into the physical characteristics of rotating BS matter can be obtained from the energy-momentum 
tensor [251] 
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(2.33) 



in which the shear tensions Tg^ and Tr , as well as all main tensions are different from each other, i.e. Tr^ ^ Tg ^ T^"^, 
and non- vanishing. Compared to a rotating ideal fLuid of a neutron star with T^"^ — {p + p)u^u^ —pg^'^, there are no 
shear tensions at all and the main stresses generally obey T^*" = Tg'' ^ T^'^ . 

A rotating BS, i.e. with large self-interaction constant A » 47rm^/Mp[, was investigated by Ryan [238]. For the 
construction of the solutions, the fact that the derivatives of the scalar field vanish (due to the scaling with A) yields 
an ideal fluid of BS matter as in the non-rotating case [58]. From the solutions, multipole moments like mass, angular 
momentum, and mass quadrupole moment are determined. For fixed angular momentum J, the maximal mass was 
calculated by increasing the maximal scalar field value until the numerical computation became unstable (Fig. 2 in 
[238]); under this condition, the maximal mass is about 0.58 \/AMpj/m^ in our units which are different by a factor 
^/A■K to the ones in [238]; further details in Section IIIF. 

More recently, in flat spacetime soliton rings are constructed [11] in a similar way from a complex scalar. They are 
metastable due to the same combination of conserved Noether charge and angular momentum. 



It has been demonstrated that BSs can actually form from a primordial bosonic "cloud" by Jeans instability 
[161,92,116,28,193,288,152,185,135,235]. (The primordial formation of non-gravitating non-topological soliton solu- 
tions was studied by Frieman et al. [101,102].) That the quantum aspects of this ground state instability of the 
complex scalar fleld acts like a classical Jeans instability is shown in [117]. The formation of a spherically symmetric 
mini-BS as a Bose-Einstein condensation is considered in [30] . One of the first steps in the formation process of CBSs 
is investigated in [71,72] by calculating the dynamics of charged bosonic clouds; cf. Section IIIF. 

In a numerical simulation, Seidel and Suen [269] have shown for a spherically symmetric mini-BS, that there exists 
a dissipationless relaxation process calling gravitational cooling; cf. [266,267,18]. This mechanism is similar to the 
violent relaxation of collisionless stellar systems which settle to a centrally denser system by sending some of their 
members to larger radii. Likewise, a bosonic cloud would gravitationally cool to a BS by ejecting parts of the scalar 
matter. It seems that the influence of the viscous term in the KG equation (2.3) is not efficient enough, so that 
radiation of scalar particles is the only viable mechanism. This was demonstrated numerically by starting with a 
spherically symmetric configuration with Mjnitiai > AfKaup, i-c which is more massive than the Kaup limit, and then 
monitoring its time evolution during which the star oscillates and scalar matter is radiated. Actually, such oscillating 
and pulsating BS branches were predicted earlier in our stability analysis via catastrophe theory [171,248]; cf. Section 
II F. Oscillating soliton stars were constructed as well by using real scalar fields which are periodic in time [268]; 
cf. Section IV. Without imposing spherical symmetry, the emission of gravitational waves have to be included. In 
Ref. [74] , it is argued that the purely radial collapse considered in [269] is too artifically. A primordial bosonic cloud 
could have had too large angular momentum so that the minimal densities required for BS formation via gravitational 
cooling are not reached. 

The situation of Bosc gas clouds, partly even before Bose-Einstein condensation, was investigated in [137] and [77]. 
Details of [137] are exhibited in [152] already, so that we concentrate on the scenario of a non-isothermal Bose gas 
cloud [77] in Section VI in which BS formation in different perfect fluid approximations are outlined. 

The influence of a cosmological phase transition on a BS with a field quantised real scalar is considered in Section 



The temperature dependence of the soliton star of Friedberg et al. [99] is investigated in [279] . In the calculation for 
the case without gravitational coupling, i.e. for a non-topological soliton (NTS), a critical temperature of = 0.7875cro 
is derived and, in the following, it is assumed that the influence of gravitation does not change Tc very much. 
Furthermore, the zero temperature estimates of [99] can be applied as well up to a temperature of T = 0.7fTo. 
Without gravity, the mass of the NTS vanishes at Tc, prohibiting their existence for T > T^, the same holds for the 
soliton star. Below the critical temperature, radius and mass of a soliton star increase with temperature as one may 
expect since the kinetic energy of the particles grows and the potential barrier between true and false vacuum gets 
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reduced with higher temperatures. If spontaneous symmetry breaking is restored at T^, the soUton star as well as the 

NTS vanishes. 

The scenario of a mini-BS with temperature is considered in [30] . Again, a critical temperature is derived numerically 
(Tc = 0.582m): above the bosons are in a gaseous phase, at Tc there is a second order phase transition, and below 
there is a mixed soliton-gas phase. A hot mini-BS with m = 1 GeV/c^ has a radius in the order of magnitude of 
cm which is about 13-14 orders of magnitude larger than in the cold state, cf. Table IV. If the scalar field mass is 
close to the Planck mass, the condensation begins with a first order phase transition. For the case m = Mpi/10, 
Tc = 0.00607m. 

Soliton BHs at finite temperature are mentioned in Section VII. For real scalar BSs, temperature (or some kind of 
motion inside the star as in oscillating BSs) is decisive for its existence, see Section IV. 



L. Boson stars in Jordan-Brans-Dicke theory 

In the so-called scalar-tensor (ST) theories, a real scalar field replaces Newton's gravitational constant G, the 
strength of whose coupling to the metric is given by a function ^((/ijbd). In the simplest scenario, the Jordan-Brans- 
Dicke (JBD) theory [38], zu is a. constant. GR is attained in the limit 1/zu 0. To ensure that the weak-field 
limit of this theory agrees with current observations, w must exceed 500 at 95% confidence limit [312] from solar 
system timing experiments, i.e. experiments taking place in the current cosmic time. This limit is both stronger and 
less model-dependent than bounds from nucleosynthesis [2]. Scalar tensor theories have regained popularity through 
inflationary scenarios based upon them [175], and because a JBD model with zu = —1 is the low-energy limit of 
superstring theory [96,46,189]. 

In this part, we exhibit the BS theories which contain beside the complex scalar field of the BS matter, the real 
scalar field of a JBD type theory. We point out the differences in the Lagrangians, but hope also to lighten the links 
between the papers. The results lead to observable distinctions between the models. In the BS matter part, the 
authors use in general the potential ?7csw so that a scalar field self-interaction can be taken into accoimt. 

The Lagrangian for our system of ST gravity coupled to a self-interacting, complex scalar field in the (physical) 
Jordan frame is 



yjBD 



R 



(2.34) 



The gravitational scalar is (^jbd and to(^jbd) is the Jordan frame coupling of ^jbd to the matter. In more general 
theories, the real scalar 0jbd possesses even a potential V . The complex scalar $ has mass m and is self-interacting 
through the potential term ?7csw = t^CSw(I^P)- 

There is an alternative representation of Lagrangian (2.34) in the so-called Einstein frame. The transition to this 
frame is effected by the conformal transformation 



Qfjiv — ^ ^^^^ 9ixu 



(2.35) 



where 



(2.36) 



and a((pE) is the Wagoner transformation [303] from 0jbd to the gravitational scalar ^p■e, in the Einstein frame. The 
relationship between 'cu{(l)}^-a) and a{ip^) is obtained by requiring 



da 



, 9(/3e J 2tj7 -|- 3 

Using the potential V{ip^) = e'^'^'^'^'^'y [(/)jbd(¥'e)]) we find the Lagrangian in the Einstein frame 



(2.37) 



2k 



+ 



_SMg2a(yE) 



(2.38) 



Mathematically, the transition from the Jordan to the Einstein frame is a conformal change of metric, cf. [202], 
and a field redefinition or "renormalisation" of the scalar field, the so-called Wagoner transformation. The question 
which "frame" describes the true, physical metric that measures the distance between spacetime points, is a subtle 
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one which depends on the couphng to matter, see the careful analysis of Brans [39]. However, the Einstein frame is 
sometimes used for calculations of BSs and in certain situations the results are close enough to GR. In the following, 
we will indicate which frame has been implemented by the authors. 

Both the coupling functions ^(^jbd) in the Jordan frame and the a((^E) in the Einstein frame are a priori unknown. 
There are, however, some theoretical reasons to motivate explicit forms. Furthermore, once a particular form for the 
coupling functions is taken, there are experimental constraints that can bo imposed. 

In Table V, we refer to publications which consider different gravitational theories and the frame in which the 
calculations are done; in some papers, more than one version is investigated. 

TABLE V. Publications on complex scalar field BS models in scalar-tensor theory calculated 
either in the Jordan or Einstein frame. 



Scalar-tensor theory 


Jordan frame 


Einstein frame 


Non-minimal-gravity coupling 
JBD, vj = const, V = 
JBD, ru = const, F 7^ 
ST , = ^(^jbd), a = a{ipE) 


[29] 

[284,119,293,294,309-311] 
[292,293,310,311] 


[59,16,15] 

[315,93] 

[59] 



BS with non-minimal gravity coupling, 0jbd = 1 -I- 2k5|$|^, w = V = 0: In 1987, van der Bij and Gleiser [29] 
considered a special JBD case, for which the gravitational scalar is actually a function of the BS complex scalar field, 
i.e. the BS matter changes its own gravitational field. The effective gravitational constant in the coupling constant 
l/(167r(jeff) ■= <?!'jbd/(2k) = l/(167rG) -|-^|$|^ can become infinite if the free parameter ^ is negative. The following 
critical masses and particle numbers arise 



^ ^ +00 , M„it ~ 0.73 V C Mii/m , N„n « 0.88 V ^ M^i/m^ 
^^-00, Merit « 0.66^ M|i/m, N„n ^ 0.72^\ M^Jm^ 



(2.39) 
(2.40) 



BS in JBD theory and coupling between (^jbd and ^, tu = Kp /A, V = ^ to^^jbd: In [284], Tao and Xue 

applied a JBD theory with positive constant vj including a free parameter /. It allowed for a space-dependent boson 
mass, determined by the JBD gravitational scalar, the corresponding term in the complex scalar field part reads 
't^^';^jbd1^|^- After redefinition, w = Gp /A, calculations are done for = 0, 1,2 • 10^ and a self-interaction constant 
A = 0, 2, 20, 60 (redefined in our units, see Appendix). For p = 1, the critical values are 



M, 



crit 



0.17\/AM|i/m , TVcrit « 0.18\/AM|i/m^ . 



(2.41) 



BS in JBD theory, zu = &,V = Q: Gundersen and Jensen [119] used a pure JBD theory with a potential C/csw and 
07 = 6 is chosen for the early universe. The critical mass is found to be 



Merit « 0.213\/AM|i^„ek/m . 



(2.42) 



Thus, the deviation of the critical mass is just a few percent from the GR limit value 0.218. 

BS in ST theory in Jordan frame, txj = tj7(0jbd), V = 0: In 1997, Torres [292] selected three different functional 
forms for the JBD coupling function, namely 



1. 2zu{(j)jBT>) + 3 : 



2Bi 



2. 2k7((/)jbd) + 3 = tUo^jBD , 

3. zu{x) = O.lx , lOx , ln{x 



|1 - 0jbd/</'jbd,o|" 

exp(0.01a;) , 



(2.43) 

(2.44) 
(2.45) 



where Bi,a,n are all positive constants, wo a constant, and 0jbd,o is the asymptotic value of the inverse of Newton's 
constant (e.g., the present value). Choice 1. had already been investigated in a cosmological context [24], for case 
2. an analytical solution is known [23] and a constraint from nucleosynthesis was determined [291]. The intention 
for explicit functions in situation 3. was to understand how the behaviour of va can be managed within the radius 
of the star. Calculations have been done in the Jordan frame for several combinations of Bi = 2, 5, 8 together with 
a = 0.5, 1.0, 1.5, 2.0 and A = 0, 10, 50, 100, 150, 200 or roo = 2 and n = 3, respectively. Again the order of magnitude of 
the BS mass varies only about a few percent when compared with GR. The influence of different external asymptotic 
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values of the gravitational scalar i^jbd on the BS mass has been investigated as well with the result that the mass 
decreases with increasing time. 

BS in JBD/ST theory in Einstein frame, a = a((^E), V = Q: Comer and Shinkai [59] studied BSs in Einstein frame 
calculations using two forms for a{ip-s), one being the Brans-Dicke coupling 

«(</'e) = ,^ , (2.46) 

with the constant Wjbb = 600, and the quadratic coupling 

«(</'e) = \{^l- Veoo) , (2.47) 

which is a particular form considered by Damour and Nordtvedt [66] except for the additive constant; in this cosmo- 
logical model there is a constraint of 77 > 3/8. The term (/Jeoo represents the asymptotic value of the gravitational 
scalar field. Tests in the solar system provide the constraint Wjw > 500; observations of binary pulsars require 77 > —5 
[69]. BS calculations are performed for r] = 0.38 which is close to the limit 3/8. In the Jordan frame, the coupling 
function is 

2xu{(I>jbd) + 3 = \ , . (2.48) 

In both cases, solution sequences are produced for A = 0, 10, 100 and for zero, one, and two nodes in the complex 
scalar field solution; hence, we find here the first excited states for BSs in scalar tensor theories. Again, there is no 
significant hut still observable difference to GR in the critical mass values. Interesting is that there is a minimal value 
of the central ipe{^) for which BSs can exist. Even more, a non-miiquencss is found, i.e. two different mass values 
exist for one '/?e(0), and an intersection point at which these two masses are identically. Furthermore, the stability of 
such BSs are discussed by using catastrophe theory following [171,172,248,173,174]. Constructions of BSs in the early 
universe proved that the BS mass decreases if we go back in time as Torres [292] noted before. In the early universe, 
only BS states with positive binding energy did exist, meaning that no stable BS could have formed. This result was 
later doubted by Whinnett [309] as we are going to mention below. 

BS in JBD theory in Einstein frame, w = 60, 600; 1^ = 0: Numerical simulations for the dynamical evolution of 
spherically symmetric BSs have been carried out by Balakrishna and Shinkai [16,15]. They also investigated the BH 
formation of a perturbed equilibrium configuration on an unstable branch. A perturbed unstable BS can also migrate 
to a stable solution by emitting scalar waves; it is shown how an excited BS state passes into the ground state (similar 
to the GR case calculations some years earlier [267,18]). Furthermore, the formation of a stable BS from a Gaussian 
scalar field packet is demonstrated in JBD theory, in much the same way as in GR. 

BS in JBD/ST theory in the Jordan frame, vu — g7(0jbd), V = 0, gravitational memory effect: The main focus 
in [293] was the gravitational memory effect of BSs both in JBD theory with w = 400 and in ST theory with 
2ix7((^jbd) + 3 10/-\/|l — (0jbd/</'jbd,o) both in the Jordan frame; this ST theory is a special choice of case 1. in 
Eq. (2.43) [292]. 

When the gravitational coupling is evolving, this has important implications for astrophysical objects [22], because 
it means that the asymptotic boundary condition for (^jbd is a function of epoch affecting the structure of the BS: (a) 
The star can adjust its structure in a quasi-stationary manner. Two sub cases can be distinguished: The adjustment 
time tadj is smaller than the evolution time tu of the Universe; then the BS will always be in equilibrium with the 
external development. If fadj > ^U) the BS will be (quasi-)static, hence, "conserving" the value of G at formation time. 
This effect is called gravitational memory. Actually, a BS which is not completely decoupled from the cosmological 
expansion will develop slowly, (b) The BS interior produces a feedback on the asymptotic gravitational constant. If a 
high density of BSs formed, they might be able to reduce the gravitational interaction strength in quite a significant 
region around themselves. 

BS in JBD theory in Jordan frame, vu = ~1, V = 0, tensor mass: In [309], Whinnett analyses three different 
proposed mass functions for BS systems for the low-energy limit of superstring theory where n7 = — 1 [96,46,189]. 
These are the Schwarzschild mass corresponding to the ADM mass in the Jordan frame, the Keplerian mass in the 
Jordan frame and the Keplerian mass in the Einstein frame which he calls tensor mass; cf. with the earlier investigation 
in [176]. In his Figure 1, the results for the mass definitions are compared with the results for the particle number in 
the Jordan frame. The mass definitions differ significantly leading to contrary physical interpretations for stability. 
The Keplerian mass in the Jordan frame yields positive binding energy for all values of central density, suggesting 
that every solution is generically unstable. The Schwarzschild mass in the Jordan frame instead leads to negative 
binding energies for every value of central density, suggesting that every solution is potentially stable. However, the 
tensor mass peaks at the same location as the rest mass, an important property in the application of catastrophe 
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theory to analyse stability properties [171]. Therefore, Whinnett adopts the tensor mass as the real m,ass of the star. 
We stress the interesting situation that whereas the (Keplerian) mass is calculated in the Einstein frame, the particle 
number (or rest mass) is evaluated in the Jordan frame. Finally, we remark that, in this paper, the dilaton field if is 
used, i.e. (/)jbd = exp(-(^). 

A comparison of the critical masses shows that for the choice w = 6 of Gunderson and Jensen, the mass values are 
about 7% different from the results of the tensor mass [119]. Already for w = 500 the difference between the masses 
is less than 10~^, i.e. negligible. 

Whinnett [309] proved also that the tensor mass and the particle number peak at the same value of the central 
density, which, in GR, is known due to Jetzer [151]. In the same year, Yazadjiev [315] gave a similar analytical proof 
for more general ST theories with a term -F(^jbd)-R, cc7(0jbd) = -f^(0jBD)</'jBD, and F(^jbd) ^ where F,H,V are 
general functions. 

BS in JBD theory in Jordan frame, zu = 400, V = 0: More physical investigations for such BSs at different cosmic 
epochs are carried out in [294] . It results that BSs can be stable at any time of cosmic history and that equilibrium 
stars are denser in the past applying the tensor mass. It is shown that the radius corresponding to the critical BS mass 
remains roughly the same during cosmological evolution. Several new physical features are displayed: the mass-radius 
relation, the behaviour of the difference between the central and asymptotic (/)jbd value, the dependence of the binding 
energy on the coupling constant zd. 

Charged BS in JBD/ST theory in the Jordan frame, V = 0: In 1999, Whinnett and Torres [310] investigated the 
influence of an additional ^7(1) charge for such BSs with constant coupling function w and a special choice of power- 
law ST theory. Solutions are calculated for the following parameter combinations (always (/'jbd,oo = 

1): For A = 

w = —1, 1, 500, 16(/)jBj,/3, for A — » 00 K7 = 1. It is found that there is a maximal charge per mass ratio, q = eMpi/m, 
above which weak field BS solutions are not stable. For q = 0, the mass evolution for four different stable BSs is shown 
for tu = — 1 and N =constant if (ijBD.oo changes from 1 to 6 (Fig. 8 therein); for different constant charges q the same 
is done in their Fig. 9. For BSs with constant central density it is demonstrated that the BS mass increases under 
certain conditions for the BS solutions. Additionally, the effect of spontaneous scalarization is investigated [311]. 

BS in JBD/ST theory in Jordan frame, V = 0: In [311], Whinnett considered the appearance of spontaneous 
scalarization. This phenomenon can also develop in neutron stars as shown in [67] in the Einstein frame. It means 
that for a star below its critical mass, the scalar tensor field is nearly constant throughout the star but, for more 
massive stars, (pE has a large spatial variation. Besides BSs with w = 1, 10, 500, 3300, the function a{(j)E) = —k4>E is 
investigated in the Einstein frame which is equivalent to the Jordan frame coupling 



2t^(<^jBD) + 3 = , . (2.49) 

2k\og(j)jsD 



The third choice is 



where fc = 1,3 and ^ = 1,2 are the constants considered and (/)jbd > 1 in the third choice. For all three cases, the 
complex scalar field self-interaction constant is set to 0; additionally, for the JBD case and the coupling function (2.49), 
A ^ oo. The calculations are carried out in the Jordan frame. Whinnett concludes that there is no spontaneous 
scalarization if A = 0, but it appears in the situation with a strong self-interaction, and so, ultimately depending on 
the compactness of the star. 

BS in JBD theory in Einstein frame, w = Q^V — V{4>-e)'- Fiziev et al. [93] examined the influence of different poten- 
tials V giving a mass to the dilaton field 0e, cf. [88]. The chosen potentials are: (a) V — (3m^jj/2)[l — exp(20E/\/3)]^ 
which corresponds to V = (3m^ii/2)(0jBD — 1)^ in Jordan frame, (b) V = 2TO^;[(/)g, (c) V = 2m^;[ sin^(0E), (d) 
V = 4mjj[[l — cos(0e)], (e) V = 2m\:^^[l — exp(— 0g)]. The parameter 7 = mm/m is introduced and all calculations 
exhibited in the paper are performed for potential (a) and A = 20 in our conventions; it is stated that potentials 
(b)-(e) does not change the results noticeable. It is discovered that for large values of 7 ~ 10 the GR solutions are 
foimd. An increasing 7 means an increasing dilaton mass with respect to the boson mass leading to a smaller dilaton 
range and, hence, suppressing the dilaton field. Additionally, the stability analysis for the BS states is performed 
using catastrophe theory. 



M. Boson halo 



In 1994-96, two different models for fitting rotation curve data of spiral galaxies were proposed using BS-like objects 
[275,157,250,251,178,254-256,259]. In [275], solutions with n = 5 and n = 6 of the Newtonian BS system (2.28), (2.29) 
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were applied to fit the ripples of spiral galaxy NGC2998 without the contribution of visible matter. The influence of 
visible matter was included in [157] for NGC2998 and NGC801 by assuming some functional matter distributions for 
the bulge and the disk. In [178], the rotation velocities for a BS with A = and A — 300 for n — % are shown without 
fitting rotation curve data. In all three papers, just the Newtonian rotation velocity formula t;^(r) = M(r)/r was 
applied, i.e. the non-negligible influence of the radial pressure of the anisotropic BS matter was ignored; the correct 
formula [255] reads 

^2(^) = ^ + ^p,(r)r2e^('')+'^('-) . (2.51) 

In Refs. [250,251], Newtonian solutions for the massless case of the Lagrangian (2.1) were investigated. Later these 
solutions were designated as hoson halos [255]. It was shown [250] that the boson halo mass increases linearly with 
radius so that the rotation velocity is constant with distance; both oscillate slightly, e.g. the rotation velocity around 
a constant value as 

(2.52) 

In [251], massless self-interacting complex scalar fields show a similar behaviour. However, for too strong gravitational 
fields, independent on the self-interaction constant, a singularity forms at the centre; cf. [199] where also the case of 
a massless conformally coupled complex scalar field with U = i?]$]^/6 leads to a linearly increasing mass. 

A more detailed investigation of the boson halos included the influence of the visible HI gas and star matter 
[254-256]. The spiral and dwarf galaxies fltted in [255] belong to a sample of 11 galaxies of different Hubble types and 
absolute magnitudes which fulfil several strong requirements [26]. All rotation curve data are measured in the 21-cm 
line of neutral hydrogen, so that the gas distribution extends far beyond the optical disc (at least 8 scale lengths) and 
the necessity of dark matter is becoming obvious. Isolation of galaxies are another constraint so that perturbative 
effects of nearby situated galaxies are negligible. For dwarf galaxies, the dark matter density near the centre is almost 
constant, i.e. dark matter has a core in these galaxies. The Newtonian solutions of [255] for a massless complex scalar 
field reveal a constant density near the centre as well, and so, good fits for dwarf galaxies could be obtained. In 
all fits for dwarf galaxies, it is recognisable that the dark matter halo dominates the luminous parts of the galaxies. 
Especially the maximum (and the following drop) of the DD0154 data can be matched perfectly by the boson halo; 
i.e. decreasing rotation curve data, near the end of observational resolution, can be explained by using a dominating 
dark matter component. The sample of rotation curve data of spiral galaxies are fitted as well as the universal rotation 
curves of [230]. Furthermore, oscillating behaviour within the optical rotation curves of some galaxies can be found, 
e.g. in NGC2998, which can be fitted by a dominating boson halo part. The general relativistic configurations for the 
boson halo reveal very large gravitational redshifts (calculated up to a value of 20). In [17], it is shown that boson 
halos are stable. 

Recently, it has been shown that BS solutions with a soliton type potential, U = m^]<I>]^(l — A6l<I>l^), can be used 
as well to fit rotation curve data of spiral and dwarf galaxies [216,217]. Further investigations on boson halos with 
a massive complex scalar field derive a scalar field mass of about 10~^^ eV/c^, if a self-interaction is included the 
scalar field weighs about 1 cV/c^, cf. [8,9]. Cosmological consequences are considered in [135,240]. A prediction on 
core radii from a formation scenario is given in [235]. 

In a different approach [75,76], the dark matter part consists of a gas of massive bosonic particles, and a total of 
36 galaxies were fitted including also a contribution for the baryonic matter. The best fits for all galaxies is achieved 
for a mass of about 60-70 eV/c^ of the bosonic particles. 

N. Classical Klein-Gordon particle 

Conventionally, a BS turns out to be an object of a huge number A'' of particles, but there are also some attempts 
to construct a BS with just one particle [91,236]. Feinblum and McKinley [91], shortly before the work of Kaup, 
investigated the mini-BS case with two constraints, namely the normalisation of both the particle number N = lo 
and of the mass M = Lo'^/m. This leads to exact one ground state particle solution, possessing a vanishing binding 
energy by construction. The authors had numerical problems to find a normalised ground state; this may be due 
to the chosen mass value, rn^ = 0.1, which seems not to be the appropriate value of the single ground state, as 
Kaup has pointed out [159]. We may add that, in fiat spacetime, the KG equation does not permit a one particle 
interpretation, as is well-known. In geometrodynamics, there are different approaches to describe elementary particles 
[308,200,201,205,207]. 
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In 1994, Rosen [236] constructed particle- type solutions for the linear KG equation which he called Klein- Gordon 
particles. The first constraint is that the total mass is equal to the particle mass, i.e. M = m. The second constraint 
is rather ad-hoc, namely m = oj. Furthermore, it is assumed that these particles have a clear surface where the scalar 
field is zero and the vacuum Schwarzschild metric starts to be present, similarly as for neutron stars. The mass of 
the particles are in the order of the Planck mass. In a follow-up paper [236] , Rosen constructs Proca particles with 
the same conditions as above. Due to an oscillating electric potential, necessarily also an oscillating magnetic vector 
potential is needed for this construction. In Ref. [199], it was noted that the radial pressure for the KG and the Proca 
particle does not vanish at the surface, thus these solutions are unstable. 



O. £)-dimensional boson stars 



Whereas in four dimensions only numerical BS solutions are known, a complex scalar with large self-interaction 
constant A » 1 allows analytical solutions in (2+1) dimensions. Spherical symmetry reduces in (2+1) dimensions to 
circular symmetry. Such circular BSs in a spacctime with negative cosmological constant exist and can be stable as 
long as the absolute value of the negative cosmological constant is lower than some critical one, otherwise all BSs have 
positive binding energy and are not energetically favourable. This result was derived for rotating BSs as well, and for 
non-rotating boson- fermion stars [241-243] . Furthermore, the rotating BS shows a vacuum hole at its centre, similar 
to the (3+1) dimensional counterpart, cf. Section III. An analytical solution for a Newtonian (3+l)-dimensional BS 
in the limit A 3> 1 is determined in [178], cf. Section II G. 

Recently, mini-BS solutions in three, four, or five dimensions with an asymptotic spacetime with negative cosmo- 
logical constant has been constructed in Ref. [10] . The effect of the negative cosmological constant on the BS interior 
is an additional attractive interaction. Therefore, the mini-BS mass decreases with decreasing negative cosmological 
constant. Whereas a numerical investigation of the radial perturbations shows that stable mini-BSs exist, from the 
Figures, it follows that all BSs have positive binding energy and, hence, are potentially unstable. 



P. Gauge boson solitons 

There exists a huge amount of papers on self-gravitating SU (n)-solutions, i.e. for Einstein- Yang- Mills theories, and 
combinations with other matter fields [25,41,42,302]. These globally regular solutions can be regarded as gauge BSs, 
cf. [249], except that they are unstable [278,43]. Abelian fields occur in the gravitational- electromagnetic entities 
{geons) of Wheeler [307] ; cf. [6] . If the gravitational infiuence is weakened so much that even Newtonian gravity can 
be neglected, non-topological or topological solitons are obtained [97,179], cf. [4], or for cylindrical symmetry, — in a 
way — thread 'BSs' can be found [221,222]. 



Q. Axidilatons from effective string models 



Commonly, in four-dimensional effective string theories [68] , the tensor field g^i, of gravity is accompanied by several 
scalar fields. Besides the familiar dilaton (p, the 'universal' axion a, a pseudo-scalar potential for the Kalb-Ramond 
(KR) three-form H := e^/^-^ *da arises. Through spontaneous compactification from ten dimensions onto an isotropic 

six torus of radius e^ , a further modulus field /? emerges. 



In the string frame g^^i, = e "^^^f g^i,, the effective string Lagrangian reads 
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(2.53) 



cf. Eq. (11) of Dereli et al. [79]. Via conformal change [202] of the metric = exp((/3/\/2K)gp,y, we can go 

over to the Einstein frame. Thereby, the kinetic dilaton term changes sign and allows us to formally combine [272] 
the axion and the dilaton into a single complex scalar field 



^:=a-\-if^e-'*'/f^ , 

the axidilaton. For this complex scalar field the Lagrangian reads 



R_ 
2k 



(2.54) 



(2.55) 
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The emerging self-interaction has, in lowest order, a mass term 

U{ip, a) := e'^'^/f^U{ip, a) ~ | $ |2 , (2.56) 

which exhibits residual C/(l) symmetry, resembling symmetry restoration. Thus, the U{1) sector of this effective 
string model is invariant under the global Noether symmetry $ — > e*"$ which allows us to establish again the global 
stability [171,172] of the star. 

All the occurring masses m and decay constants / of the bosonic particles are related via m^f^ ~ mafa — TOtt/tt — 
10^^ (cY)'^ to those of the pion. Whereas the mass of the dilaton Lp is conventionally related [70] to the supcrsymmetry 
breaking scale msusY by m^p ~ 10^''(msusY/TcV/c^) eV/c^, this is not the only possibility. In the context of string 
cosmology [105], massless axions arc able to seed the observed anisotropy of the Cosmic Microwave Background 
(CMB). The same holds for the KR axion, provided it lies in an ultra-low mass window. Besides the full dilaton 
interaction for the bosons as mentioned above, we could consider a very light dilaton if which is stabilised [81,82] 
through the axion. 

Instead of combining the dilaton with the axion, one could also take recourse to the other moduli field and consider, 
e.g., the complex Kdhler form field ^ := ip -\- i\/vl(3 akin to T-duality [83], without including the axion. 

Are MACHOs axidilaton stars? In Rcf. [213,214], we have proposed that the weakly interacting dilatons and 
axions, after a partial Bose- Einstein condensation to an axidilaton star (ADS), account also for some fraction of 
the MACHOs (Missive Compact Halo Objects) which are dark matter candidates, almost certainly detected by 
gravitational microlensing [281], though still with statistical ambiguities [118]. Since the total mass of an ADS is 
not well constrained, one could phcnomcnologically turn the argument around: By identifying the MACHOs of 
gravitational mass of about 0.6 Mq (as given by microlensing observations) with the critical mass of an ADS, we 
are essentially "weighing", via MxaupZ-^crit — m, the constituent mass to be of the order nia — 10~^° eV/c^. It is 
encouraging to note that such an ultra-low mass is perfectly compatible with the constraints on the mass range of 
the KR axion seeding the large-scale CMB anisotropy, cf. the recent results of Gasperini and Veneziano [105] within 
low-energy string cosmology. 

If ADSs during their evolution would accumulate a larger mass, they can start to oscillate [171] and thereby get 
rid of some excess mass due to 'gravitational cooling' [269]. Repeated accretion leads ultimately to a BH in the 
upper range of the MACHO mass. These axion induced BHs [271,104] do not carry scalar hair, which could serve to 
distinguish them from primordial BHs, but could have some remnant of P-violation or even "axion hair" [84] , if a is 
interpreted as a supcrpotcntial for the KR axial torsion. In the pre-big bang scenarios, the dilaton would produce a 
detectable gravitational wave background [61]. 

Therefore, if such string-inspired scalar fields would exist in Nature, axions could not only solve the non-baryonic 
dark matter problem [298] , but their gravitationally confined axidilaton stars would also represent the MACHOs which 
presumably account for about 20% of the halo dark matter in our Galaxy. 



R. Quantised complex scalar field boson star 

Field quantisation of a complex scalar has been studied in [117] in the context of BS formation. In Ref. [132], a BS 
with a field quantised complex scalar has been introduced. 

As is well-known, the complex scalar $ can be regarded as two real fields, namely (f>i = + $*)/\/2 and (1)2 = 
i($ — $*)/v^. Then, field quantisation results in the operator expansion 

<I>J = E^n?J^F^"Wi^J"'(^, ^)e-^-f^?* + 6li?P>(r)[yl«l(^, ^)]*e+^-f^^ , (2.57) 

nla 

where j = 1,2. For a bound state, i?p''"(r) are radial distributions as generalisations of the wave fimction of the 

hydrogen atom, vj^'^O, ip) = (l/-\/47r)P;''''(cos6') e"'!"!''' the spherical harmonics given in terms of the normalised 
Legendre polynomials, and |a| < ? are the quantum numbers of azimuthal a.nd an^M^ar momentum. The non- vanishing 
commutation relations for the bosonic creation and annihilation operators are 

&Ml']=Snn'Su,Saa'. (2.58) 

In general, there arc two number operators N^^'> := b'^^^'^b^^K For the ground state of a cold BS, if A''i = N2, we have 
|A^) = \N, 0, 0, . . .) := ^io(/|0) = T\i ^io(/|0), where |0) is a vacuum state in the curved spacetime 'background'; 
cf. Section IV A. 
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In Rcf. [132] an attempt was made to motivate symmetry breaking in a model with |<I>|^ self-interaction and self- 
gravity. Three different BS models lead to three different effective coupling constants. The BS with a not quantised 
complex scalar [58] has simply A/2. The BS model with a field quantised real scalar produces 3A/4; and the BS 
model with a field quantised complex scalar has 5 A/8. All three models results in the limit A 1 to small changes 
in the critical mass Merit ~ K^fkM'^Jm, where K = 0.22, 0.225, and 0.223, respectively. 

III. HOW TO DETECT COMPLEX SCALAR FIELD BOSON STARS 

In this section, we speculate on observational consequences for the spherically symmetric BS model given by (2.1), 
if not mentioned differently. In several cases, we shall recognise the important influence of the scalar field potential. 
Under our assumption that the complex scalar field has no interaction with baryonic matter other than gravitationally, 
the BS is transparent, i.e. baryonic matter can accumulate up to the centre of the BS or photons and particles can 
move through the BS interior, respectively. 

To detect a BS depends on the physical situation; there is a corresponding wave-band where the BS might be very 
luminous. In the following, we shall discuss several of such scenarios. 

A. Rotation curves 

It might be possible to detect a BS in X-rays. Imagine a massive BS, say of IQ^Mq, for which it is likely that an 
accretion disk forms and since its exterior solution is asymptotically Schwarzschild, nearly beyond its effective radius 
the BS looks similar to an Active Galatic A^ucleus (AGN) where usually a black hole (BH) at the centre is assumed. In 
X-rays [252,258,264], one can probe close to the Schwarzschild radius or, as we propose, even inside the BS; Iwasawa 
et al. [139] claimed that by using data from the satellite ASCA they have probed to within 1.5 Schwarzschild radii; 
of. [177,197,320]. Because this is inside the static limit for a non-rotating BH, they conclude that a Kerr geometry is 
required. A BS configuration could still be an alternative explanation, giving a non-singular solution where emission 
can occur even from the centre. 

In order to differentiate between a BS and a BH, it is interesting to calculate the tangential velocities of baryonic 
matter that may rotate around the centre of a BS. For the static spherically symmetric metric (2.7), geodesies of a 
collisionless circular orbit obey 

vl = \ry'e^ = le^'ie' - 1) + f P.r^e^+^ - ^ + f P.r^e^+'^ , (3.1) 

which reduces for a weak gravitational field into the Newtonian form j^^^^ = M{r)/r if Pr = 0. Rotation curves for 

the cases A = 0, 10 and 300 were calculated for a critical mass BS in [252,258]. The curves increase from zero at the 
centre up to a maximum followed by a Keplerian decrease. The maximal rotation velocities reaches still inside the 
BS more than one-third of the velocity of light; similar results attain for the potentials Ucg, Usg, Ul (Fig. 8 of [262]). 
For a BS in Newtonian approximation, the maximal circular velocity can be just about 2000 km/s and this at more 
than 50 times of its Schwarzschild radius, but still inside the BS. One can conclude that large rotation velocities 
are not necessarily signatures of BHs. Baryonic matter rotating with the maximal velocity of about c/3 possesses an 
impressive kinetic energy of up to 6% of its rest mass. If one supposes that each year a mass of IMq transfers this 
amount of kinetic energy into radiation, a BS would have a luminosity of lO^^erg/ s. 

B. Gravitational redshift 

If baryonic matter inside a BS emits or absorbs radiation within its gravitational potential, the spectral feature is red- 
shifted [252,258,264] . The gravitational redshift Zg within a static background is given by 1 + Zg := ^ ^^{R^^i) j ^v{Ri^t) ^ 
see e.g. Ref. [219], where emitter and receiver are located at K\-a\, and i?cxtj respectively. The maximal possible redshift 
for a given configuration is obtained if the emitter is exactly at the centre i?int = 0. The receiver is always practically 
at infinity, hence exp [^'(iiext)] = 1- For all other redshifts in between, the gravitational redshift function is defined as 

l + ^g(a;) = exp(-^) . (3.2) 
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Tabic 1 in [252] gives the results for the critical mass BS for different A. The redshift lies between 0.0657 for A = —20, 
0.4565 for mini-BS with A = 0, and 0.6873 for the limit case A 3> 1. With increasing self-interaction coupling constant 
A also the maximal redshift grows. Unstable BSs presumably may have arbitrary gravitational redshift. In Ref. [262], 
for the cosh-Gordon potential Ucg, the maximal redshift is found to be 0.46, for the C/(l)-Liouville potential U^, it is 
-Zmax = 0.49. For a neutron star, Markov [196] derived Zmax = 0.49. In general, observed redshift values would consist 
of a combination of cosmological and gravitational redshifts given by (1 + ^) = (1 + Zc){l + Zg) due to the addition 
formula of relativistic velocities. 



C. Gravitational lensing 

In this section, the results of gravitational lensing of a transparent spherically symmetric mini-BS are shown [65,263]. 
In contrast to the last sections, it is assumed that the BS interior is empty of baryonic matter, so that deflected photons 
can travel freely through the BS. The deflection angle [304] is then given by 

oo 

J rv r^e — 

ro 

where b is the impact parameter and rg denotes the closest distance between a light ray and the centre of the BS. 
The lens equation for small deflection angles can be expressed [218] as 

sin(t?-/3) = ^sina, (3.4) 

where D\s and Dqs are the distances from the lens (deflector) to the source and from the observer to the source, 
respectively. The true angular position of the source is denoted by /?, whereas i? stands for the image positions. One 
usually defines the reduced deflection angle to he a = d — (3 = sm~^ (Dis sina/Dos) • However, equation (3.4) relies 
on substitution of the distance from the source to the point of minimal approach by the distance from the lens to the 
source Dig, see [218]. For large deflection angles the distance D\s cannot be considered a constant but it is a function 
of the deflection angle so that the form of the lens equation changes [65] into 



sin a = — — cos v cos 



arsin ( sm{'d — a) 



[tana? + tan(d-t?)] . (3.5) 



The reduced deflection angle can be kept deflned as a = i? — /3. 

Numerical computations of the reduced deflection angle for different BS potentials Uk, UqGi Usg, Ui, were performed 
by assuming that the BS lens is half-way between the observer and the source, i.e. Dis/Dos = 1/2 (for details, see 
[65,262]). Observable differences are fomid depending on the choice of the self-interaction. In the case of a simple 
mass term Uk, the largest possible value of a is 23.03 degrees with an image at about i? = n x 2.88 arcsec where 
n = n{Doi,uj) is the distance factor which is a function of the distance from the observer to the lens and the scalar 
field frequency the inverse of which is associated with the BS radius. For non-relativistic BS approximations, smaller 
angles will occur. The angle of the image position can have very different orders of magnitude, depending on n. For 
example, n = I fixes to be measured in arcsec. Under the assumption that the BS mass is 10^*^ Mq, the distance 
Do\ is about 100 pc. If the distance factor is n = 10~^, then ^ is measured in milli-arcsec and the BS-lens is at about 
100 kpc from the observer. 
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FIG. 2. Reduced deflection angle for a mini-BS of critical mass. The angle for the image positions is given in units of 
arcsec times the distance function n [65]. 

In Ref. [65], it is summarised that the BS has all qualitative features of a non-singular spherically symmetric 
transparent lens, cf. Fig. 2, the lens curve for a mini-BS of critical mass. Three images can be observed, two of them 
being inside the Einstein radius and one outside. An Einstein ring with infinite tangential magnification, the so-called 
tangential critical curve, is found, and also a radial critical curve for which two internal images merge. The appearance 
of the radial critical curve distinguishes BSs from other extended and non- transparent lenses [246]. For a BH or a 
neutron star, the radial critical curve does not exist because it is inside the event horizon or the star, respectively. 
Two bright images near the centre of the BS and the third image at some large distance from the centre are found. 
If an extended transparent source like the BS is considered, one can imagine how to detect the BS. In such a case, 
two images, one radially and the other tangentially elongated and both close to each other, might be observed. By 
looking along the line defined by these two images, the third one can be detected at a large distance. 

For the Cosh-Gordon potential Ucg [262], the maximal reduced defiection angle is 23.229 degrees and for ?7l an 
even larger deviation, namely 24.391 degrees, with an image at about the same place is found. Differences among all 
four mentioned potentials are up to 4%. 

Discussions of lensing effects of different objects can also be found [239,169]. 

D. Gravitational microlensing and MACHOS 

A statistical analysis of the galactic halo via microlensing [226,281] suggests that MACHOs (massive compact 
halo objects) accovrnt for a significant part (> 20%) of the total halo mass of our galaxy. Their most likely mass 
range seems to be in the range 0.3 - 0.8 Mq, with an average mass of 0.5 Mq, cf. [155,281]. If the bulge is more 
massive than the standard halo model assumes, the average MACHO mass [155] will be somewhat lower at ~ 0.1 
Mq. However, there are some astrophysical difficulties with this result of estimated ~ 0.5 Mq for the lenses. These 
cannot be hydrogen-burning stars in the halo since such objects are limited to less than 3% of the halo mass by 
deep star counts [114]. Modifying the halo model to slow down the lens velocities can reduce the implied lens mass 
somewhat, but probably not below the sub-stellar limit 0.08 Mq. Old white dwarfs have about the right mass and 
can evade the direct-detection constraints, but it is difficult to form them with high efficiency, and there may be 
problems with overproduction of metals and overproduction of light at high redshifts from the luminous stars which 
were the progenitors of the white dwarfs [49]. Primordial BHs are a viable possibility, though the coincidence has to 
be explained to have them in a stellar mass range. 

Due to these difficulties of getting MACHOs in the inferred mass range without violating other constraints, there 
was the suggestion that BSs could be the explanation [213]. The explicit physical values for a mini-BS has been 
exhibited in Fig. 3 with a scalar field mass of 10~^° eV /(?. 
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FIG. 3. Left: Mass M and particle number N (or rest mass mN at infinity) of a mini-BS depending on the central density 
p. Right: Mass-radius dependence of a mini-BS. 

At the left-hand side of Fig. 3, the dependence of the mass M and the particle number N (or rest mass) on the 
central density p is shown. Stable non-rotating BSs exist at the lower central densities below the maximum mass. 
The critical values arc Afcrit = 0.846 Mq and mA'crit = 0.873 Mq for a central density of pc ~ 9.1 x pnuci, where 
Pnuci = 2.8 X 10^^ kg/m'^ is the average density of nuclei. Since non-interacting bosons are very "soft" (due to 
Bose-Einstein condensation), BSs are extremely dense objects with a critical density higher than comparable ones of 
neutron or strange stars [110]. The figure on the right-hand side gives the mass depending on the radius. For the 
mass-radius diagram, 99.9% of the total mass was chosen as effective radius. This ensures that the exponentially 
decreasing exosphere of the BS has almost no influence on the asymptotic Schwarzschild spacetime. 

Stable BSs have radii larger than the minimum at 20.5 km with a mass of 0.846 Mq. The choice for the boson mass 
m yielded that the total mass of these relativistic BSs is just in the observed range of 0.3 to 0.8 Mq for MACHOs. 

As was the case in the discussion of axidilaton stars (Section HQ), one could turn this argument around: By 
identifying the most massive MACHOs with known gravitational mass of about 0.8 Mq as BS, one is essentially 
'weighing', via MKaup/A^crit — m, the boson mass to m 10"-^'^ eV/c^. 

E. Cerenkov radiation from boson stars 

As in Section III C, we assume that the BS interior is empty, but now electrically charged particles move fast through 
the gravitational BS field. In general, if a particle moves through a medium with a constant velocity greater than 
the velocity of light in that medium then Cerenkov radiation occurs. Because of the superluminal particle motion, a 
shock wave is created and this yields a loss of energy. 

Accelerated charged particles emit electromagnetic radiation, but acceleration only by a gravitational field does not 
produce Cerenkov radiation, in agreement with the equivalence principle. There, the acceleration induced by gravity 
disappears in the local inertial frame (further details in the references of [47]). Despite that, the gravitational BS field 
can act as a medium with an effective refractive index for light if particles move with constant velocity through it 
[120]. It is found 

n^(fco) = 1 - , (3.6) 

ftp 

where ko is the frequency of the emitted photon 7 and R'^^ is the sum on the spatial indices of the Ricci tensor. 
Cerenkov radiation is now kinematically allowed if i?'^ < 0, i.e. n^{ko) > 1. 

Following [120], it has been investigated in [47] whether BS models allow Cerenkov radiation. It was found that 
stable mini-BSs and stable BSs cannot generate Cerenkov radiation, whereas unstable BSs can. However, non- 
topological soliton stars do allow this effect where the position and the strength depend on the central scalar field 
value and the potential parameter $o- Then, Cerenkov radiation appears in spherical shells and only there, but this 
can be everywhere inside the star. 

Again, we recognise how important the infiuence of the scalar field potential on detection can be. 
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F. Gravitational waves 



In the early stage of BS formation, highly excited configurations are expected to exist in which the quantum numbers 
n, I and a of the gravitational atom, i.e. the number n — 1 of nodes, the angular momentum and the azimuthal angular 
dependence e"*"*^ are non-zero; cf. [148]. For energetic reasons, these excited BSs are presumably only meta-stable, 
all initially higher modes during the BS formation have eventually to decay into the ground state n = 1, Z = (with 
a = for non-rotating and a ^ for rotating BSs) by a combined emission of scalar matter radiation and gravitational 
radiation. 

If we assume a BS to be in an excited state with Z = a = 0, we know that its mass and particle number depend 
linearly on the number of nodes n— 1, cf. [98,29]. Since these states have zero quadrupole moment, the only decay 
channel is scalar matter radiation. For the transition of an excited BS with only n > 1 into the ground state, it is 
reckoned [92] that the energy released by scalar radiation is about 

^rad ~ (n - l)M|i/m , A7V ~ (n - l){Mpi/mf , (3.7) 

where the loss of bosonic particles is given by A A''. 

Fcrrcll and Gleiser [92] considered the situation of a Newtonian BS that has the main part of particles in the Is 
state and only a small number in the 3d (0.286 parts per thousand). They estimated the amount of total energy 
released by the transition from this small number of particles in the 3d into the Is state. The lowest mode which 
has quadrupole moment and therefore can radiate gravitational waves is the 3d state with n = 3 and 1 = 2. For 
the An = 2, AZ = 2 transition, the BS will decay into the Is ground state with n = 1 and I = 0, and thereby 
the gravitational radiation process preserves the particle number N. Even though only a relative small number of 
particles populates the 3d state, the radiated energy is quite large [92], i.e., E'l-ad = 2.9 x 10^^ (GeV/mc^) Ws. The 
typical frequency of this gravitational radiation is t'g = 2.03 x 10^^ (mc^/GeV) Hz with the power of 1.3 x 10'^^ W 
but much less than the output of a supernova explosion of about 10^^ W. Highly relativistic BSs with radius near the 
Schwarzschild one (M/ R 1/G) could have a maximal power of 10^^ W, but this is just a theoretical upper bound 
on the gravitational radiation of a self- interacting system. In the scenario of [92], it is concluded that the excited 
state will decay extremely fast to the ground state with the exponential decay time of about 1.1 x 10~^^ (mc^/GeV) 
s. Thus, the final phase of the BS formation would terminate in an outburst of gravitational radiation despite the 
smallness of the object, cf. Table IV. 

In [72], the collapse scenario of parts of a charged boson cloud with uj = m, i.e., just above the eigenfrcquencies of 
excited CBSs, is investigated. In comparison with the results above, gravitational waves are created in a shorter time 
interval (10~^^ s) and with much less power (4 x 10^^ W). 

The issue of non-radial pulsations of a BS was mathematically formulated [166], whereas quasinormal modes are 
explicitly calculated in [316]. The interest in non-radial pulsations arises from the fact that they transport gravitational 
waves away from the BS. For linearised perturbations around an equilibrium BS, odd and even parity modes can be 
distinguished. It was shown that the odd parity modes do not couple to gravitational waves as well as for the stellar 
pulsation of a perfect fluid. This mode type manifests the propagation of the gravitational wave through the star, 
it does not change the star's density or pressure distributions. With a different analytical approach, the even-parity 
quadrupole oscillations {1 = 2 modes) which do couple to gravitational waves were investigated numerically in [316]. 
These quasinormal modes describe oscillations with emitted gravitational and scalar matter radiation (but no incoming 
radiation). In contradistinction to relativistic fluid stars, the imaginary parts of the frequencies are large, i.e. the / 
and p modes imply very short damping time scales. It is also found that the amplitude of the metric perturbation 
is comparable to the one of the scalar field, again different to the w modes of an ordinary star. No weakly damped 
modes are present within the BS. One restriction of the numerical calculations in [316] is connected with the definition 
of the BS surface. Due to numerical reasons, a clear surface had to be defined, but a BS has only an exponentially 
decreasing exosphere; this has a crucial influence on the values of the eigenfrequencies; an introduction to pulsating 
relativistic stars can be found in [167]. 

Future gravitational-wave detectors may have also different ways of searching for BSs. If a compact object of, say, 
a solar mass is observed to be spiralling into a central one with a much larger mass, one would be able to test the 
larger object. From the emitted gravitational waves, the values of the lowest few multipole moments of the central 
object can be extracted, like mass M, angular momentum J, mass quadrupole moment M2. For a BH as central 
object, the no-hair theorem states that all its multipole moments arc uniquely determined by the two lowest ones, 
e.g. M2 = —J^/M. In [238] M, J, M2 and the spin octopole moment 5*3 are calculated for rotating BSs in distinction 
to a BH. For example, if one observes the combinations J/M^ = 0.01, -M2M/J2 = 24, and SsM'^/J^ = 19, then the 
massive central object would be very likely a rotating BS. If instead, one flnds — M2M/J^ = 1 and S^M'^/J^ = 1. 
then a BH is the likely source, and if one determines — M2M/J^ = 24 and S^M^ / = 4, a different configuration 
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(BS with different kind of self-interaction, e.g.) may be detected. Further numerical combinations for BSs result from 

the Figures in [238]. 

For a rotating BS, co-rotating particles orbiting with circular geodesies in the equatorial plane were investigated in 
[238]. The particle energy depending on its gravitational- wave frequency for one BS solution is shown there in Figure 
6. If particles could travel freely inside the BS, it would, in principle, be possible to map out the interior of the BS 
via the emitted gravitational waves. 

If a BS experiences some kind of perturbation, we may be able to measure characteristic frequencies. In [123], 
the mode frequencies for the fundamental mode and the first harmonic mode of radial perturbations of mini-BSs are 
exhibited in Tables II and III (cf. Figs. 7 and 8 as well); cf. Fig. 6 in [267] where similar frequencies were discovered. 
Further details of Ref. [123] are mentioned in Section II F. 

The output of gravitational waves was considered in [108] for identical BS binary systems, assuming that they 
are located in galactic halos. However, the gravitational radiation background contains data from ordinary binary 
systems as well. The estimation for critical BSs gives a contribution to the gravitational radiation background where 
the amplitude is /i ~ 3 x 10^^'"' (TeV/mc^)'^/^ and the frequency v ~ 25 (TeV/mc2)5/4 Hz. This result is valid if 
the BS binary density is about the critical density and they dominate over ordinary binaries if m < 10^ TeV/c^. 
For mini-BS binaries with scalar field mass m = 10~^ eV/c^, the typical amplitude is /i ~ 10"^'' and the frequency 
V 4.5 Hz, again for a mini-BS binary density of about the critical density. 

G. Boson star in the galactic centre? 

Several galaxy centres contain 10^ up to 10^ solar masses. The widely accepted explanation is a supermassive BH 
or a rich cluster of objects (stars, BHs, e.g.). Already in [252,258], a so-called mini-BS or BS has been proposed as 
an alternative explanation, cf. Sections III A and IIIB. 

In [295], a supermassive mini-BS as model for the Galactic centre was compared with the observational data 
available for our Galaxy where, at least, the existence of a star cluster can be excluded [85,86]; cf. [106]. The central 
object has a mass of about 2.6 x IO^Mq. The more recent observations of [106] probe the gravitational potential at 
a radius larger than 4 x 10^ Schwarzschild radii of such a BH. The mentioned data consist of the movement of stars 
around the Galactic centre and, in [295], they were explicitly fitted by a solution of a mini-BS with m, ^ 10^^^ cV /c? 
plus a stellar cluster. The same result can be received by using a BS solution with m ^ GeV/c^ (A = 1) and a 
non-topological soliton star solution (<I>o ~ m) with m ~ 10** GeV/c^. The data can be fitted by a central BH as well 
and, so far, no distinction can be obtained between a BH and a BS at the centre. The difference starts at a radius 
more than three orders of magnitude below the innermost data point, in case of a general relativistic mini-BS [295]. 
A mini-BS in the Newtonian approximation has less mass, but if we reduce the constituent mass m in order to obtain 
the same central mass, the radius becomes about two orders of magnitude larger (cf. Fig. 5 and Table 2 in [252]). 
Observations are just one to two order of magnitudes away from finding differences in the BH-BS comparison. Let us 
also mention that an extended neutrino ball can explain the data [297] so far. 

A non-rotating BS may have also the problem of accreting too much matter and so forming a BH eventually. It 
is obvious that inspiralling interstellar gas and stars will collide with each other, glue together within the BS, and 
so eventually form a BH. If instead a small BH spirals inwards and stays there, the final state of a non-rotating BS 
will be without doubts a BH. The investigation of the effective potentials for massless and massive particles of a 
spherically symmetric mini-BS can give some understanding how infalling matter behaves inside the BS [295]. All 
particles with orbital angular momentum and unbound orbit do not move through the BS centre, hence, leave the 
star; only if the orbital angular momentum is zero the particles will move exactly through the centre. In this ideal 
situation of non-interacting matter, every matter is removed from the BS interior and no BH can form. But. as 
mentioned above, presumably matter will collide and glue together so that a BH cannot be avoided (i.e., the particles 
have bound orbits). However, in case of a rotating BS, there is the possibility that matter could be guided along 
geodesies on the rotation axis forming jets; this has still to be shown. For a BS in the limit A ^ 1, circular orbits are 
stable until about 5/3 times the radius of the doughnut hole [238]. 

On the basis of an effective potential, it is shown [296] that, for a spherically symmetric mini-BS, circular orbits 
exist for every value of the radial coordinate and are stable, even inside the BS. Hence, accretion could follow a series 
of stable orbits and finally end up at the centre. Further investigations are needed in order to clarify whether a BH 
forms. 

What happens if the BS under external perturbations (gravitating scalar matter) is forced to leave the equilibrium 
ground state changing mass and particle number in the configuration? The numerical simulations of [267,269,18,15] 
show that stable mini-BSs under finite non-infinitesimal perturbations start to oscillate, emit scalar field radiation, 
and settle down into a new mini-BS configuration with less mass and a larger radius. If the amount of scalar accretion 
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matter does not exceed a critical value during some time, the BS is not jeopardised to form a BH. 

The emissivity properties of a geometrically thin, optically thick, steady accretion disc with constant accretion rate 
around a supermassive critical mini-BS are demonstrated in [296] and compared with a BH of same mass. Optical 
thickness means that the local effective temperature must be sufficient to radiate away the local energy production. 
The disc calculations are only valid where the gravitational potential is ~ GM/r, i.e. outside of the investigated 
mini-BS. One result is that, for the inner disc parts, a BS produces more power per unit area and a hotter disc 
than the comparable BH; for this critical mini-BS, there is a temperature difference of about two thousand degrees 
in the inner parts of the accretion disc. For a BS with A = 100, there is about the same maximal temperature but 
the position of that maximum is closer to the BS centre. The calculation of the final emission spectrum reveals a 
difference between the mini-BS and the BH model in the far ultraviolet regime > 10^^ Hz (Fig. 3 in [296]); this disc 
model shows no production of X-ray or gamma-ray energy. The Eddington luminosity, at the balance of the inward 
force of gravity with the outward pressure of radiation, is not exceeded within the mini-BS interior. Thus, accretion 
could continue inwards a mini-BS. A further result [296] is that the disc model does not describe the Galactic centre 
Sgr A* neither for a BH nor a BS due to the accretion disc properties itself. It produces several order of magnitudes 
too much accretion luminosity. In contrast to observations is also the peak of the standard accretion disc broad band 
spectrum in the infrared. 

Additionally, BSs could disrupt stars [295]. In general, there is a distance to a central object where the extension of 
a normal star cannot be neglected and where tidal forces on the star become important. The corresponding so-called 
tidal radius is defined where M/r' is equal to the mean internal energy of the passing star. Estimations show that 
BHs should have masses below IO^Mq so that they provoke disruption outside the event horizon. In general, such an 
event takes place once in about 10** yrs. Hence, whereas a too massive BH disrupts stars just inside its event horizon, 
a BS produces this effect visible for observations (inside or outside of its effective radius). 

Observations in the x-ray band within the central regions of active galaxies have shown an interesting feature, a 
Doppler and gravitationally redshifted iron Ka line [283,89]. It was claimed [139] that the data can be understood as 
being emitted from within 1.5 Schwarzschild radii near a rotating BH. In [252,258], we speculated that a BS geometry 
may be able to explain the data as well. Recently in [190], a detailed numerical analysis of a geometrically thin, 
optically thick accretion disc around a supermassive critical mini-BS simulated this situation for different inclination 
angles. It is concluded that the observed line profile of the galaxy MCG-6-30-15 cannot be fitted by a spherically 
symmetric mini-BS or BH for the chosen accretion disc model. For a possible identification of such a non-rotating 
mini-BS, the characteristic spectral feature would be the detection of an intense double peak for face-on galaxies or, 
for larger inclination angles, of several peaks of the long drawn-out Ka line. 

H. Boson star inside a HII cloud? 

BSs could exist in two more scenarios as it was pointed out in [252,258]: 

Within host galaxies of quasars, bright Hll regions arc observed [14]. By imagining that neutral Hi gas clouds could 
gain kinetic energy in the gravitational potentials of BSs, an excited hydrogen gas cloud, i.e. Hll, results. 

The second situation is that a light BS is completely contained within a nuclear burning star. In this case, both 
stellar striictures would be interacting mutually, especially changing the model parameters of the conventional star. 
This resembles the so-called Thorne-Zytkow object, i.e., a neutron star core inside a supergiant [286]. If instead the 
BS matter has some kind of interaction with high energetic photons, it might be possible that the BS dissociates 
and the remaining scalar particles circulate freely within the conventional star. A so-called cosmion represents such 
a scenario and tries to explain the solar neutrino problem [90]. 

I. Gravitational memory and the evolution of boson stars 

The evolution of conventional nuclear burning stars follows the Hertzsprung-Russell diagram of the luminosity as 
a function of the temperature. During its life time, such a star has different stages which are expressed by a line in 
this diagram. There is the possibility that a BS could go through a continually evolution as well, if the gravitational 
attraction changes during the evolution of the Universe. 

In JBD or scalar tensor theory, a real scalar field regulates the strength of the gravitational force, cf. Section II L. 
Suppose that at the time of BS formation, the gravitational strength was different from the one prevailing today; 
following this theory, G was larger in the past. The real scalar field is actually a radial function so that the gravitational 
strength inside the JBD-BS is determined by a mutual gravitational interaction of the complex and the real scalar 
field matter. These investigations are outlined in Section II L. From the observational side, it may be important that, 
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inside the BS, still the strength of the gravitational constant at formation time could be effective [293]. According to 
this gravitational memory effect, the BS is book-keeping the evolution of G. Because the JBD field is a radial function 
and changes its value at infinity, there should be a repercussion on the BS. But if this change is much slower than 
the cosmological evolution of G, the star is practically static. Due to this effect, physical properties like the radius 
depend on formation time and could be different even for the same total mass. 

If instead the JBD field adapts quickly inside the JBD-BS as changes occur at infinity, the BS evolves as well. 
Figure 1 in [294] shows the mass against the value of the real scalar field at infinity for constant particle number. It is 
derived that the JBD-BS mass increases with increasing (/)jbd(oc) = 1/G, i.e. with time. One can understand that an 
increasing JBD field pumps energy into the BS and increases the mass. Thereby, the evolution of a BS in JBD theory 
demonstrated in this Figure 1 resembles the evolution in the Hertzsprung- Russell diagram of conventional stars. The 
physical properties as mass, radius, and stability is influenced by this repercussion [294], such that, for example, the 
mass changes about 2 %. 



IV. REAL SCALAR FIELD BOSON STARS 



In contrast to the complex case, real scalar fields do not possess a conserved charge or particle number. Moreover, 
a static configuration $(r) = P{r), i.e. (2.6) for uj = 0, would, in fiat spacetime, be unstable due to Derrick's theo- 
rem [80], whereas a superposition of positive and negative frequency states (/)i(r, = (P(r)e~*'^* -|- P(r)e*'^') /V2= 

\/2P{r) cos{ujt) is time-dependent and implies time-dependence of the energy-momentum tensor and the metric, as 
well. These configurations arc called oscillating soliton stars [268,269] or 'oscillatons' [300,301] and are, in general, 
unstable. Actually, our global stability analysis via catastrophe theory [171,172] has revealed already earlier the 
existence of oscillating BSs, as we will call them in the context of this review. Q-stars are a class of non-topological 
solitons which are already discussed in detail in [152]; cf. [57,191,270,12,280]. Further interesting configurations are 
the so-called pulsons [32,33] or oscillons [134], respectively; these are localised, time-dependent, unstable, spherically 
symmetric solutions to a nonlinear KG equation. 

The temperature T is the decisive parameter characterising the thermodynamic equilibrium of such a boson gas of 
real scalar particles [187]. The particle density in momentum space is given by 

= eME^/kBT) - 1 ' ^^-^^ 

where = yj^+rr? is the energy of the particle with momentum p and mass m. In the limit of vanishing 
temperature, one obtains n$ = 0. Hence, a BS consisting of real scalars cannot exist for T = and a different method 
has to ensure a conserved particle number. 

For a massless real scalar [55], analytical solutions of the Einstein-KG system are derived first by Buchdahl [45], 
cf. Wyman [314] and Baekler et al. [13]. Wyman showed that all static spherically symmetric solutions are contained 
in the two-parameter family of Buchdahl. The common characteristic of all such solutions is that they are singular 
or topological non-trivial [163,164]. In [244,245], the ADM mass of such solutions were calculated because the metric 
is asymptotically flat and the scalar field vanishes at infinity. Since the latter has a logarithmic singularity at the 
centre, for some solutions, the metric is singular as well and these are typical examples for naked singularities. 

By including a mass term and a quartic self-interaction, numerical solutions of the Einstein-KG system can be 
determined [153]. Again the scalar field possesses a singularity at the origin; alternatively, one finds regular solu- 
tions at the centre, but the singularity is shifted to a finite distance. The dynamical stability investigation used a 
variational perturbation method estimating upper bounds for an eigenvalue. The result showed that all such naked 
singularities are unstable [153,156]. However, a direct numerical determination of the ground state eigenvalue indicates 
perturbatively stable solutions [56] . 

In the numerical onc-parameter-solutions of Ref. [153], the 'bare' mass m is the decisive parameter, whereas A is not 
important. In the case of complex scalars, two parameters occur: besides m, we have the frequency w, and globally 
regular solutions exist. Thus, the parameter w in the phase factor is actually responsible for the regularity of the 
solutions. This is not surprising because this parameter is proportional to the conserved particle number A^. This 
consideration motivated the construction of the globally regular one-parameter-solutions with non-vanishing uj and 
m = in [250,251,254-256]; cf. boson halos in Section II M. 
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A. Quantised real scaleir fields 



An assembly of real scalar particles experiencing just a mass term was already considered in 1966 by Bonazzola and 
Pacini [34]. In this approach, the scalar field is second quantised. Later on, sohitions for this model of "systems of 
self-gravitating" bosons were considered by Rufiini and Bonazzola [237] similarly to the Hartree-Fock atom. Because 
a real scalar has no anti-particle states, the corresponding KG field in a spherically symmetric spacetime metric can 
be merely decomposed into positive and negative frequency field operators 

$ = $+ + $-, (4.2) 

where 

nla 

= 5^6„iai?r(r-)[V;'"'(^, <p)]*e+-'''* (4.3) 

nla 

are again generalisations of the wave function of the hydrogen atom for a bound state. Here Rf{r) arc radial distri- 
butions, Yi"'^{0, ip) — (l/-\/47r)P;'"'(cos^?) e^*!"!''' the spherical harmonics given in terms of the normalised Legendre 
polynomials, and |a| < I are the quantum numbers of azimuthal and angular momentum. (The nodeless radial field 
solution RQ{r) is earlier denoted by P{r).) The non- vanishing 6osonic commutation relations are 

[Kla, b\'l'a'] = ^nn'5w5aa' ■ (4.4) 

For the ground state of a cold configuration, 

N 

\N) = \N,0,0,...):=l[^\oo\0) (4-5) 
1 

is chosen [237], where |0) is a vacuum state in the curved spacetime 'background'. The BS energy- momentum tensor 
Tpi,($) becomes now an operator. Hence, for the right-hand side of the Einstein equation (2.2), the vacuum expectation 
value := {N\ : T^i, : \N) is calculated for the ground state where : T^^, : denotes normal ordering of the operator 

products. 

An excited state of positive energy can be defined by \N, n, I, a) := $+10). Such a 'gravitational atom' [92] represents 

a coherent quantum state, which nevertheless can have macroscopic size and large mass. The gravitational field is 
self-generated via the mean value of the energy-momentum tensor, but remains completely classical, whereas the real 
scalar is treated to some extent as operator. 

As we noted in the beginning of this Section, there is no conserved particle number for real scalars, and so one 
needs to introduce the normalisation condition 

oo 

27r J r^ef^-'^)/^ [^^^R^i(^r) + cOn'i'RnH'ir)] dr = Snn'Sw (4.6) 



in order to stabilise the system. Similarly as for complex scalars, the constraint is imposed that the solutions are 
asymptotically flat at spatial inflnity. 

Due to the fact that the mean value of the energy-momentum tensor is employed, the same system of difi'erential 
equations (2.2), (2.3) arises as for the complex case. A later investigation [40] rediscovered these results. 

An interesting issue is how a BS reacts to a second-order phase transition in the early Universe. In [132], such a 
spontaneous symmetry breaking has been investigated by assiiming a negative quadratic term in the scalar potential. 
After the phase transition, the 'bare' mass m is increased and the self-interaction strength is reduced by A ^ A/2, 
resulting in a smaller total mass M of the BS and, thus, preserving its existence at all. 

In [128,129,188], the BS of a field quantised real scalar, the model of Ruffini and Bonazzola [237], has been extended 
by including a neutron star in perfect fluid description. The details are discussed in the 1992 reviews. 

The situation of a field quantised complex scalar is exhibited in Section II R. 
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B. Dilaton star 



The most general scale-invariant Lagrangian density for the real dilaton field ip coupled to gravity adopts the form 

-Cdii = vl^e^^//- I i^i? + \g>'''{d^<p){d^^) - jjy/f- 1 , (4.7) 

where A is a coupling constant and is the 'decay constant' of the dilaton. The dilaton mass itself is given by rridii = 
2\/Xf^. For this non-minimal coupling, scale invariance is realized by the conformal change g^^ exp(— 2a)(ji^;^ 
and the shift Lp ^ ip -\- af^p in the dilaton field, where a is some dimensionless parameter. If = -y/S/STrmpi, the 
theory is invariant under the larger symmetry group of conformal transformations, induced by a space-time dependent 
parameter a = a{x). For a scale-invariant interaction, the potential V{(p) = e'^'^^^fXf^/A is the only allowed type. For 
any other dilaton potential, the theory is still of ,JBD type in terms of = fipe'^^^^ , but not anymore scale invariant. 
In the Einstein frame, the theory consists just of a masslcss real scalar plus a cosmological constant A = coupled 
to gravity. 

Since the Lagrangian is scale invariant, Noether's theorem provides the conserved dilaton charge 

Qdii = (1 + 6/K/2) J ^\e^^li^d'^ d^x . (4.8) 

It vanishes if p) is time-independent as in the case of a static space-time of a non-rotating star. Therefore, in order 
to construct static dilaton stars, Gradwohl and Kalbermann [115] followed the semi-classical approach of Ruffini and 
Bonazzola [237] and Breit et al. [40] for real scalars. They regularise the energy-momentum operator by means 
of the subtraction := — (0|Pju|0) of its vacuum expectation value. In curved spacetime, this Wick or normal 
ordering of the operator products, should be performed by the point-splitting method, in order to avoid ambiguities 
due to the non- uniqueness of the vacuum state, cf. [305], p. 87. Thereby, they can effectively get rid of the cosmological 
constant. A spontaneous breaking of the scale invariance leads to a vanishing dilaton expectation value at infinity. 
As a consequence, GR is recovered asymptotically and the dilaton star carries zero dilaton charge, but has a constant 
particle number iV due to the normalisation of the dilaton field. 

The following parameter combinations have been considered: A := A(/^/a;)^ = 1, 10, 30, 100 and a = (f^/Mpi)^ = 
1/40, 3/(87r), 1, 10. The most noticeable effect on the total mass M is due to the change in the latter parameter. For 
a very light dilaton of mass mdii = 10~^^ eV/c^, Gradwohl and Kalbermann [115] found as maximal values 

McHt = 7y/X^MQ , iJerit = 40\/a"^ km , (4.9) 
where oj is the ground state frequency. 

C. Axion star 

In general, axions have no effective cooling mechanism; they remain in the form of an extended gas cloud after 
separation from the Hubble flow. Decay processes may help in cooling down an axion cloud. In a cosmological 
context, there is a temperature dependence of the axion decay modes [232]. The stimulated decay process aa ajj, 
together with aa — > 77 and aa — > uD have been examined for their importance in the cooling of an axion cluster 
[103]. (The on-shell contribution aa — > aa followed by a — > 77 can be ignored because this does not yield an energy 
loss.) As later shown in [287], the stimulated processes can only help in forming axion stars, by including the effect 
of large phase-space density; this process relaxes the axion clouds efficiently to form axion stars. On the other hand, 
free axions have a lifetime which is longer than the age of the universe so that the axion stars can be re-powered by 
accretion. Axion stars could be detectable as coherent cosmic masers [287]. In the same publication, an axion star 
for a real scalar field a with potential m^a^ -|- Aa^ is investigated resulting in a critical mass in the order of magnitude 
of Mq. Hence, for a real scalar field, we obtain the same result as Colpi et al. [58] in the case of complex scalars. 
For this potential, the axion star may radiate as a radio source. Additionally, the relaxation time of a gravitationally 
bounded cloud of axions as well as the luminosity of the axion star is estimated in [288] . 

The formation of axion stars at the QCD epoch was considered in [133]. The initial isocurvature perturbations 
yield small axion miniclusters of mass about 10^"* kg with a diameter of the order of 10^^ m. In later evolution, the 
miniclusters undergo standard hierarchical gravitational clustering. For the minicluster density of about 10^° kg/m^, 
axion annihilation aa 77 is not important, but in axion stars with densities of about 10^^ kg/m^, it makes the 
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configuration unstable [269,74]. More recently, the formation of axion halos is discussed in [135]; cf. the constraints 
on scalar parameters in [235]. 

Numerical studies [168] including nonlinear effects in the evolution of inhomogeneities in the axion field around the 
QCD epoch can lead to very dense axion clusters with densities high enough so that the stimulated decay process 
supports to lead to Bose-Einstein condensation, and eventually, to axion stars. Femto- and picolensing by axion 
miniclusters is investigated in [169]. 

A different way how axion stars may send signals are exhibited in [140-146] . There, the oscillating soliton star model 
[268] is taken as example for axion stars. Since axions couple to the electromagnetic field, axion stars could dissipate 
their energies in magnetised conducting media surrounding white dwarfs or neutron stars. Oscillating configurations 
may generate monochromatic radiation with an energy equal to the axion mass. If axion stars collide with a white 
dwarf, the latter is heated up and a detectable amount of thermal radiation can be emitted during the collision [141]. 
Even stronger effects can be expected if the axion star collides with a neutron star. If the axion star dissipates its 
whole energy in a single outburst, this could be a possible mechanism for the observed gamma ray bursts [142]. For 
example, a collision with a neutron star with a strong magnetic field of about 10^ Tesla can produce gamma rays 
up to 10^^ eV. More details of such collisions are given in [143]. A jet of baryons and leptons with Lorentz factors 
larger than 100 are generated. If the collision does not lead to the immediate destruction of the axion star and 
consists actually of several repeated collisions between both stars, this could be the origin of time-dependent complex 
properties of gamma ray bursts. A prediction for this model is the emission of radio waves with the frequency given 
by the axion mass as a byproduct to the gamma ray bursts. In Ref. [94,95], gamma ray bursts are explained by the 
relativistic detonation of electro-dilaton stars. 



V. BOSON-FERMION STARS 

In this part, we discuss the possibility of a BS located inside a fermion star or vice versa, depending which is the 

dominating component. As an example of a fermion star, a neutron star is normally considered as indicated by the 
equation of state. So far, there are no field quantised complex scalars involved. Since this topic was already discussed 
in earlier reviews we shall be rather brief here, concentrating on new publications, and show which type of scalar fields 
are applied within the combined star. 

A. Fermion soliton star with real scalar field 

A combination of a BS with a fermion star has first been examined in the context of soliton stars [182]. There, 
a semi-classical BS with the potential J7($) = m^$^(l — $/$o)^/2 for a real scalar secures the existence of non- 
topological solitons, even in the absence of gravity. Hence, only a conserved fermionic particle number Nf exists; that 
of real bosons is not well defined. Furthermore, the real scalar interacts with the fcainion field -0 through —f^j^j^, 
where / is a Yukawa type coupling constant. For simplicity, one adopts for the fermion mass rrif = /$o so that the 
fermion has zero effective mass inside the false vacuum, i.e., inside the BS. For the fermion field, a Thomas- Fermi 
approximation is chosen. Remarkably, Nf depends as mI[^ /m^/'^<^>l on the mass, i.e., fractional as for the CBS 
(cf. Section II C), whereas the total mass Mpj/m^o is similar to that of the soliton star. Fermion soliton BHs are 
mentioned in Section VII. 

The observational properties of zero-mass fermion soliton stars are investigated in [51 53]. In Ref. [51], only the 
effects of baryons entering the star are considered if the universe temperature drops below some critical value. For the 
special choice m/ = /$o of Lee and Pang, all fermions are almost massless inside the fermion soliton star. Moreover, 
the binding energy of constituent particles decrease to zero. Hence, if nuclei enter the star, they disintegrate into 
protons and neutrons, and then both into quarks. Furthermore, the fermion soliton star contents differs in the course 
of universe evolution. Above some critical temperature, protons cannot enter (or leave), and only quark pairs are 
inside the star. Below that temperature, baryons can penetrate the star's shell and just quarks remain inside this 
so-called quark soliton star. The fermion soliton stars with only quark pairs would finally radiate away all its energy 
and cease to exist; only such stars with a net number of protons survive. These protons converts the rest energy into 
thermal energy so that for a short period of time (about 10^ years) at redshifts z larger than 4, these stars could be 
X-ray emitters. Fermion soliton stars with dominating electron part may perturb the short wavelength character of 
the cosmic microwave background radiation (CMB), causing a sharp increase in the brightness temperature. 

A quark soliton star at finite temperature including the same number of (disintegrated) protons and electrons, 
quark and e+-e~ pairs, and some so far unknown fermionic particles with ponderable mass are calculated in [52,53]; 
the result is labelled ponderable soliton star. Whereas the quarks in the pure quark soliton star are almost massless, 
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the ponderable soliton star contains fermions with ponderable mass. The lifetime of the ponderable soliton star is 
from a few weeks to several years, they have a surface temperature of about 10® K, and they radiate profusely with 
a luminosity as high as lO*'^ erg/s. Because this temperature is close to the hydrogen recombination temperature of 
~ 10^ K, these stars could be observed at redshift distances of 2; ~ 10^. This short radiation time is still longer than 
the evolutionary time of the Universe at that epoch, and so, there may be an interaction with the CMB as well. The 
expected perturbations on CMB spectrum is less than 1% and it is largely around the short- wavelength at 1 mm or 
less, in form of point radio sources; spatial inhomogeneities of arcsec scale should be observed. 

The situation of the fermion soliton star of Lee and Pang is extended in two respects in Ref. [62,63]. First, it is 
shown that the soliton star is connected to the Lee- Wick model with scalar potential 

[/($) = im2$2(i _ $/$o)2 +B 4- 3(4>/$o) («'/^o)^ • (5-1) 

For a constant i?^/"* of about 100 MeV as used in hadron spectroscopy, the fermion soliton star has a critical mass 
of about just three solar masses. Additionally, the fermion soliton star at finite temperature is investigated and by 
this the formation and evolution of the stars. Above a critical temperature around 100 MeV, the universe is filled 
by the homogeneous real scalar $ = $o- Below the critical temperature at which first-order phase transition occurs, 
the real scalar adopts its true vacuum state $ = 0; the false vacuum will only survive inside some bubbles filled with 
fermions. The bubbles contract until they are stabilised by the fermionic pressure, i.e., the fermion soliton stars have 
formed. In the star's evolution, it is assumed that the star could evaporate into hadrons either at its surface or inside 
via nucleated bubbles. Then, at present, fermion soliton stars would have a mass of about 1O~^M0. 

A real scalar obeying a linear KG equation and the above mentioned Lee-Pang interaction arc the ingredients of 
the model in [27]. The fermionic matter has zero temperature so that one should expect no real scalar matter at 
T = 0; responsible for the occurrence of the non-conserved scalar matter seems to be the Yukawa interaction. In 
the model, instead of the exact stress-energy tensor of a Dirac field as in [223], the perfect fluid approximation for 
the neutron star interior is used. Actually, Ruffini and Bonazzola [237] proved that a large number of degenerate 
fermions can bo approximated by the perfect fluid approach. Additionally, the self-interaction factor ipip is replaced 
by its expectation value following [237] . Then, the influence of the equation of state for an ideal Fermi gas for neutron 
stars is investigated, following Chandrasekhar [48], and nucleon-nucleon interactions are considered as well. It is not 
surprising that the hypothetical presence of a BS can change the cooling and the structure of the final state of a 
neutron star very effectively. 



B. Boson-fermion star with complex scalar field 

The stability of combined boson-fermion stars built from a complex scalar is investigated in [130,151]. A slowly 
rotating boson-neutron star based on Kaup's linear potential and Chandrasekhar's equation of state [48] has been 
analysed in [276]. 

In Ref. [35,36], a boson-fermion star is calculated in a generalised JBD theory with w = where the Brans-Dicke 
field has a potential and, hence, a mass; for the complex scalar, the potential Ucsw is used and solutions are derived 
for A = 0.01, 10. This work can be seen as an extension of [93]; cf. the JBD Section IIL. 



VI. PERFECT FLUID APPROXIMATION 



Several papers on BSs followed the old idea of Oppenheimer and Volkofi^ [223] and applied a fluid approximation 
for the bosonic matter, which yields the correct order of magnitude for the total mass. In this Section, the bosonic 
fluid is isotropic. 

Nishimura and Yamaguchi [220] constructed an equation of state for a cold Higgs fluid 
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where m, A arc the parameters of the Higgs Lagrangian, K, E the first and second elliptic functions, respectively, and 
u is some parameter. For several choices of w^/A, it is shown that the equation of state p{p) is an extrapolation of 
the equation of state of nuclear matter to higher density regions; a further characteristic is that the pressure vanishes 
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below some value of the central density. Substituting this equation of state into the Tolman-Oppenheimer-Volkoff 
equation leads to BSs with critical masses in the order of magnitude of a solar mass. In a way, by taking into account 
a bosonic self-interaction constant, Nishimura and Yamaguchi anticipated the result of Colpi et al. [58] who found the 
critical BS masses for the Lagrangian, i.e. anisotropic, description in the same mass region. Furthermore, these BSs 
have a finite radius, namely, the point where pressure as a function of radius vanishes. 

Takasugi and Yoshimura [282] investigated a non-self-interacting cold Bose gas with the equation of state 

p = 2,m^nlu^/^ /{2p) , p = mmu^ 1 + ^u^/s , (6.3) 

where m is the boson mass, n the number density, no some constant, and u := n/tiQ. We see that for small u, p (x p^^^ 
and for large u, p = p/3. The critical mass is found to be 0.57 Mpjm at a radius of 5.5/to which is remarkably close 
to the Kaup limit 0.633 Mpjm of the Lagrangian description. 

Three different equations of state are investigated in Ref. [299]. For masslcss bosons, the equation of state of photons 
is used p = p/3, leading to photon stars, or "geons". The radius is determined where the surface pressure equals the 
one of the cosmological background radiation. The critical mass is about 10^^ Mq and the radius is larger than the 
radius of the Universe, so that it seems that photon stars cannot be found. Alternatively, a giant MIT bag filled with 
gluon matter is considered, hence, the calculation stops where the pressure is equal to the bag constant B which is 
chosen to be 57 MeV/fm^. The equation of state can be written as p = {p — 4B)/3. Then, the critical mass of the 
gluon star is 2 Mq. The third example describes a non-rclativistic boson gas of massive particles where Bose-Einstein 
condensation is neglected, with an equation of state p = 2{p — mn)/',i. Because the pressure never vanishes in this 
model, the surface of these BSs is defined by p = B. The critical mass is about 10""^ Mq for m = 20 GeV/c^. 

The MIT bag equation of state is derived in Ref. [121] for a real scalar field with self-interaction. Thereby, the 
effective mass of the real scalar is zero, so that the corresponding BS solutions are called precarious stars. These stars 
exist only for temperatures above some value which is determined by the bag constant; this coincides with the fact 
that the real scalar has no conserved charge, cf. Section IV. 

Departing from the Newtonian spherically symmetric Poisson equation, Dehnen and Gensheimer [77] analysed a 
self-gravitating Bose gas cloud with temperature. A similar situation was already considered by Ingrosso and Ruffini 
[137] for an isothermal case, i.e., constant temperature all over the cloud; details can also be found in [152]. The 
difference from the Newtonian BSs of Section II G is that, here, an equation of state is derived that does not include 
a scalar field, but considers the Bose-Einstein statistics in its non-relativistic limit. It is found that the Bose gas 
obeys a polytropic relation p = ap'^ with some constant a. The Poisson equation reduces then to the Lanc-Emclen 
equation for polytropic index 3. In general, it is found that the temperature decreases with increasing distance from 
the centre, hence, Bose-Einstein condensation occurs in the outer regions of the star. This result is in contrast to 
the isothermal case where the condensed phase is in the centre (additionally the outer regions are Boltzmannian and 
spatially unlimited) [137]. Furthermore, these Bose gas clouds have a finite radius where the energy density, the 
pressure, and the temperature vanish. 

VII. SOLITON BLACK HOLE 

These configurations, introduced in [99], describe a non-topological soliton star sitting completely inside the 
Schwarzschild horizon which is possible due to the finiteness of such a star, cf. the radius definition in Section II D. 
That means that the radius of the star is smaller than its Schwarzschild value. The mass M in units of M|j/(m$§) 
lies between the values 0.1256 and 16/(277r) = 0.1886. For this upper bound, the particle number N vanishes. It is 
not surprising that a soliton BH is not stable, it decays very slowly by sending scalar matter towards the horizon [181]. 
Outside the horizon, the Schwarzschild solution holds; cf. the more general BS result in [228]. In [225], the physical 
properties of a soliton BH at finite temperature was investigated following the procedure for a non-topological soliton 
star in [279], cf. Section II K. Fermion soliton BHs were constructed in [182] for a real scalar field; details on fermion 
soliton stars in Section VA. 



VIII. OUTLOOK 

BSs are so far hypothetical objects; what is more it is unknown whether such localised configurations really have 
the size of a star. Due to the unknown nature of its constituents, the dimensions of BSs could be smaller than an 
atom or as large as a galactic halo. Thus, all these different scenarios have to be considered and further observational 
consequences drawn. 
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Within the last three decades, a lot of models have been labelled BSs, sometimes one and the same mathematical 
model received different physical designations. All models agree that the constituents of a BS are scalars. They differ 
in regards to whether the scalars are real or complex and whether the scalars are classical or field quantised. In 
this way, we can distinguish four basic types of BSs. Furthermore, in each of these four BS types, sub-types can be 
separated due to the form of the scalar self-interaction. In this review, we tried to clarify and sum up all different 
boson star notions. 

Further analysis is needed in order to understand these highly interesting instances of a possible fine structure in 
the energy levels of general relativistic BEC, the so-called BSs. In view of their rich and prospective structures, are 
BSs capable of explaining parts of the dark matter in the Universe? 

We hope that this review can also be understood as a starting point that the topic BS leaves the pure scientific 
area and finds its way to popular science [224]. 

We think that one cannot finish a topical review on BSs any better than by quoting from [180,98,99,182] and his 
review [184] the last two sentences of T.-D. Lee: At present, there is no experimental evidence that soliton stars [or 
BSs, the authors] exist. Nevertheless, it seems reasonable that solutions of well-tested theories, such as Einstein's GR, 
the Dirac equation, the KG equation, etc., should find their proper place in nature. 
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APPENDIX A: CONVENTION 



As usual conventions can make life easier or tougher. There are infinite possibilities to transform the BS field 

equations into dimensionlcss quantities: two are commonplace in the literature. In order to avoid recalculations, we 
would like to propose that future authors apply the following convention 

x = mr, 0=-, ^jB^a; =^^T^' '^Oz^ =^tH' ^= 'a 

m M^j Mpi 477 
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